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We theoretically consider the carrier density dependence of low-temperature electrical conduc- 
tivity in high-quality and low-disorder two-dimensional (2D) "metallic" electronic systems such as 
2D GaAs electron or hole quantum wells or doped/gated graphene. Taking into account resistive 
scattering by Coulomb disorder arising from quenched random charged impurities in the environ- 
ment, we show that the 2D conductivity o-(n) varies as ct ~ n'^'"' as a function of the 2D carrier 
density n where the exponent /3(n) is a smooth, but non-monotonic, function of density n with 
possible nontrivial extrema. In particular, the density scaling exponent /3(n) depends qualitatively 
on whether the Coulomb disorder arises primarily from remote or background charged impurities 
or short-range disorder, and can, in principle, be used to characterize the nature of the dominant 
background disorder. A specific important prediction of the theory is that for resistive scattering 
by remote charged impurities, the exponent /3 can reach a value as large as 2.7 for kpd ~ 1, where 
kp ~ -^n is the 2D Fermi wave vector and d is the separation of the remote impurities from the 2D 
layer. Such an exponent /3 (> 5/2) is surprising because unscreened Coulomb scattering by remote 
impurities gives a limiting theoretical scaling exponent of /3 = 5/2, and naively one would expect 
Pin) < 5/2 for all densities since unscreened Coulomb scattering should nominally be the situation 
bounding the resistive scattering from above. We find numerically and show theoretically that the 
maximum value of a (/3), the mobility (conductivity) exponent, for 2D semiconductor quantum 
wells is around 1.7 (2.7) for all values of d (and for both electrons and holes) with the maximum a 
occurring around kpd ~ 1. We discuss experimental scenarios for the verification of our theory. 
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I. INTRODUCTION 



It is well-known that carrier scattering by background 
quenched disorder arising from random impurities and 
defects limits the T = (i.e., low-temperature) residual 
conductivity of a metal or doped semiconductor. This 
"residual resistivity" , i.e., the extrapolated T = value of 
the electrical resistivity obtained from the measured low- 
temperature transport data, provides information about 
the extrinsic disorder in the material limiting transport 
properties. In 3D metals, this residual resistivity is not 
of much intrinsic fundamental interest except for the fact 
that very pure defect-free single crystals, where the back- 
ground disorder is presumably very low, could have ex- 
tremely low residual resistivity leading to very long elec- 
tron mean free paths. In particular, the residual resis- 
tivity of different samples (with different levels of sample 
purity, i.e., different impurity or defect content) of the 
same metal could differ by orders of magnitude, and a 
metallic resistivity uniquely defining a particular metal 
(e.g. Cu, Al, or Ag) is only meaningful at higher temper- 
atures (> 100 K) where phonon scattering dominates the 
electrical resistivity over impurity/defect/disorder scat- 
tering. At low temperatures, each metallic sample would 
have a unique resistivity reflecting its specific impurity 
content signature, and as such characterizing a metal by 
a unique resistivity is useless at low temperatures (i.e., 
each sample of the same metal would have different resis- 
tivity at T = 0). In fact, the low-temperature resistivity 



of a particular sample typically depends on the prepara- 
tion history of the sample, and annealing to room tem- 
peratures (where all samples of a particular metal do have 
the same resistivity) and then cooling down to low tem- 
peratures could substantially modify the sample resistiv- 
ity as the impurity/disorder configuration could change 
due to annealing. Theoretical statements^ about residual 
resistivity (or equivalently, low temperature resistivity) 
of 3D metals thus focus on first principles calculations 
of the quantitative aspects of the local disorder poten- 
tial arising from various defects or disorder in the metal 
(where the specific types of defect or impurity have to be 
specifically assumed), and then estimating the resistiv- 
ity arising from various postulated disorder in terms of 
resistivity due to some specified impurity or defect con- 
tent on an atomic percentage basis. Typical values of 
calculated residual resistivity for most simple metals fall 
in the ~ 0.1 — 1 /xilcm range per atomic percentage of 
impurities, and thus a metal with 99.999% (e.g., Cu) pu- 
rity could have an extremely low residual resistivity of 
10^^*^ ilcm, leading to elastic mean free paths ~ 1 cm 
although the typical phonon-limited room-temperature 
mean free paths in most metals are only 1 — 10 A with 
the room-temperature phonon-limited resistivity of most 
metals being around ^ fiflcm. 

Thus, theoretical studieg^"^ of disorder-limited residual 
resistivity of metals focus entirely on the quantitative 
modeling of various defect scattering in the metal us- 
ing detailed materials-specific band structure and Boltz- 



mann transport theories. No systematic dependence of 
the residual resistivity of metals on various metallicity 
parameters (e.g., the Wigner-Seitz radius r^ or lattice 
constant or Fermi energy) can in general be discussed 
qualitatively, and there is no metallic density scaling 
of conductivity that one can speak of since the carrier 
density cannot be tuned by an external gate voltage in 
metals (as it can be in 2D semiconductor systems or in 
graphene). In principle, there is only a modest varia- 
tion in Ts (^ 2 — 6) among different 3D metals (since 
Ts ^ n^^/-^ where n is the effective metallic free electron 
density), but the band structure variation in going from 
a metal to another would typically swamp any system- 
atic Tg-dependence of the residual resistivity. Thus, the 
only systematic trend of the residual resistivity in metals 
which is discussed in the literatureA'^ is the dependence 
of the r = impurity-limited resistivity of a particular 
metal on the atomic type (e.g. atomic number) of the 
various impurities or on the type of defects causing the 
resistive scattering. It is basically a subject of quantita- 
tive details based on serious numerical calculations. 

In the current work we are interested in discussing 
qualitatively the impurity-limited 'metallic' residual (i.e. 
T = 0) resistivity as a function of the electron density 
of the metal. In particular, we want to understand the 
density-dependence of the T — metallic resistivity as- 
suming that the metallic density can be tuned continu- 
ously keeping all other parameters (e.g., disorder, band 
structure) fixed. Of course, such a situation is practically 
impossible in 3D metals since the electron density is fixed 
in each metal and cannot be tuned at all. But, as is well- 
known, such a tunable carrier density situation is routine 
in 2D "metalhc" systems such as Si MOSFETs, gated 
semiconductor hctcrostructures and quantum wells, and 
gated graphene. In principle, a variable carrier density 
can be achieved in 3D doped semiconductor systems by 
changing the doping level (but keeping the same dopant 
element) in different samples, but this is not ideal for 
our purpose of a qualitative understanding of the den- 
sity scaling of the metallic T — saturated resistivity 
since there is always the possibility of unknown sample 
to sample variation beyond just the carrier density vari- 
ation since doping level cannot be tuned continuously in 
a single sample as can be done in 2D systems. We do, 
however, briefly consider the density scaling properties of 
3D impurity-limited saturated resistivity for the sake of 
completeness although most of the current work focuses 
on 2D semiconductor structures (and graphene) where 
the carrier density can be easily tuned in the same sam- 
ple by applying an external gate voltage, and thus the 
density-dependent conductivity is a meaningful concept 
in 2D "effective metallic" systems existing in semicon- 
ductor structures and graphene. 

Obviously, the T = conductivity (= p^^, where p is 
resistivity) would manifest different density dependence 
for different types of background disorder, i.e., different 
types of impurity-electron interaction. The main resistive 
disorder scattering in relatively pure 3D metals is due 



to defects, vacancies, and impurities which scatter pri- 
marily through the short-range (essentially) (5-function 
type scattering potential (although there are extended 
defects which could produce anisotropic longer-range dis- 
order potential). Disorder scattering by J- function type 
point scatterers gives rise to rather uninteresting den- 
sity scaling of electrical conductivity, again making 3D 
metals unsuitable for studying the density dependence of 
conductivity. Our main interest in this work is to obtain 
the density scaling of 2D T = "metallic" conductiv- 
ity arising from Coulomb disorder induced by random 
quenched charged impurity centers in the environment. 
Quenched Coulomb disorder is the dominant extrinsic 
resistive scattering mechanism in all semiconductor sys- 
tems (2D or 3D) at low temperatures since doping by im- 
purities is essential in inducing carriers in a semiconduc- 
tor (whereas metals, by definition, have free carriers at 
the Fermi level even at T = 0), leading to the inevitable 
existence of background Coulomb disorder. Thus, a key 
distinction between low-temperature transport in metals 
and semiconductors is that short-range disorder (arising 
from very strong metallic screening) dominates metallic 
transport whereas long-range Coulombic disorder (aris- 
ing from random localized charged impurities) dominates 
transport in (both 2D and 3D) doped semiconductors. 
The density dependent conductivity in doped semicon- 
ductor structures (or graphene) arises from the carrier 
screening of background Coulombic disorder which de- 
pends sensitively on the dimensionless ratio qxp/'^kp 
where qtf and kp are respectively the (density depen- 
dent) screening wave vector and Fermi wave vector of the 
system. The variation in the carrier screening properties 
as a function of the carrier density eventually leads to 
the density dependence of the resistivity. 

Of course in doped semiconductors, one must worry 
about the additional complications of disorder-induced 
Anderson localization and/or low-temperature carrier 
freezeout and/or possible percolation transition associ- 
ated with charged impurity-induced charge puddle for- 
mation. In the current work, we ignore all of these com- 
plications uncritically, focusing primarily on extremely 
high-quality modulation-doped 2D GaAs quantum well 
structures (and high-quality suspended graphene) where 
these complicating circumstances are absent down to 
very low carrier densities and very low temperatures. 
Our theoretical results presented in this paper, based on 
Drude-Boltzmann semiclassical transport theory, should 
apply to experimental systems above carrier densities 
(and temperatures) where localization (and related ef- 
fects) become operational. The neglect of carrier local- 
ization and freezeout is not a particularly significant is- 
sue for our 2D theory and results because high-quality 
2D semiconductor systems (and graphene) arc not sus- 
ceptible to these problems except perhaps at extremely 
low temperatures and carrier densities of little practi- 
cal or experimental interest. On the other hand, our 
low-temperature and low-density transport results for 3D 
doped semiconductors are given here purely for the pur- 



pose of completeness and comparison with 2D results 
since 3D doped semiconductors typically become insulat- 
ing at low carrier density (as well as low temperature) be- 
cause of localization and carrier freezeout. Our focus and 
most of our presented results, as the title of our article 
clearly indicates, are on low-tcmpcrature density depen- 
dent metallic transport properties of very high-quality 
2D systems where the concept of density scaling of con- 
ductivity is both theoretically and experimentally mean- 
ingful. 

One may wonder about the fact that the scaling theory 
of localization predicts that all 2D systems are strictly 
speaking Anderson insulators, and have, in principle, 
zero conductivity at T = for infinite samples at all 
carrier densities. (This is theoretically true even for 
graphene when disorder induced intervalley carrier scat- 
tering is taken into account.) For our purpose in the 
current work, where we are specifically interested in very 
high-mobility 2D structures with very low background 
disorder, the scaling localization is a non-issue because 
(1) the samples are finite in size and the temperature, al- 
though low, is still finite; and (2) more importantly, the 
elastic mean free path is extremely long, making the effec- 
tive 2D localization length much larger than the system 
size (or the phase breaking length, whichever is shorter). 
Thus, we are specifically considering the density depen- 
dence of the semiclassical part of the 2D resistivity in 
the situation where the semiclassical resistivity is much 
larger than the quantum (or weak localization) contribu- 
tion to the resistivity. This limit is generic in high mo- 
bility 2D GaAs semiconductor systems and in graphene, 
and therefore our work is of wide validity. Thus we are 
strictly speaking considering the density dependence of 
the zeroth order semiclassical conductivity in the situa- 
tion where the quantum contribution to the resistivity is 
negligible. 

We emphasize that the density dependence (or scal- 
ing) of electrical transport properties is rarely discussed 
in the theoretical literature mainly because (as discussed 
above) such a density dependence of electrical conduc- 
tivity is neither relevant nor interesting in 3D metals. 
Indeed, it is universally believed that the only aspect of 
electrical conductivity where studying the carrier density 
dependence is a meaningful endeavor is in the study of 
density-tuned Anderson localization as a (T = 0) quan- 
tum phase transition, exactly the physical phenomenon 
we are excluding in the current work where we focus en- 
tirely on the Drude-Boltzmann part of the conductiv- 
ity in high quality samples assuming quantum localiza- 
tion effects to be negligible. Thus, the density scaling 
of our theory is not connected with quantum critical- 
ity at all, but with the behavior of the density depen- 
dence of the background disorder arising from nontrivial 
screening properties. In contrast to the density depen- 
dence which is rarely discussed in the literature except 
in the context of metal-insulator transition, the tempera- 
ture dependence of transport properties has been exten- 
sively discussed for electronic materials (including many 



experimental and theoretical papers in 2D systems^""—) 
because the temperature dependent electrical resistiv- 
ity generically contains qualitative information about the 
underlying resistive scattering mechanism. For example, 
a linear-in-T higher temperature metallic resistivity is 
the hallmark of acoustic phonon scattering in both met- 
als and semiconductors (and in both 2D or 3D systems) 
whereas acoustic phonons typically lead to strongly sup- 
pressed high power laws (^ T"^ with ?7 sa 4 — 7 depending 
on the details) in the (2D or 3D) metallic resistivity below 
a system-dependent Bloch-Griineisen temperature Tbg- 
Our current work focuses entirely on temperatures well 
below the Bloch-Griineisen temperature (T ^ Tbg) so 
that phonons are completely ineffective in limiting the 
electrical conductivity. Optical phonons, which are often 
relevant in semiconductor transport at higher tempera- 
tures and indeed limit the room-temperature resistivity 
of 2D GaAs-based semiconductor structures^ and are ex- 
ponentially suppressed as ^-^^i^oI^bT ^i^ j^.^ tempera- 
tures (where huLO ^ 450 K in GaAs), are completely 
irrelevant for our low-temperature transport considera- 
tions. 

It is well-known that high-quality low-density 2D semi- 
conductor systems (as well as graphene) often manifest 
a strong temperature dependence in its low-temperature 
resistivity^^—, which arises from the strong temperature 
dependence of the screened Coulomb disorder^ at low 
carrier densities. Again, our theoretical work at T = 
is completely free from any temperature-induced com- 
plications in the resistivity since our interest is in un- 
derstanding the density dependence of the T = re- 
sistivity/conductivity, which can be obtained from the 
low-temperature experimental data by extrapolation to 
T = or by sitting always at a constant very low temper- 
ature (e.g., 50 mK) in obtaining the density dependence 
of the transport properties. In any case, the question we 
are interested in is perfectly well-defined as a matter of 
principle: How does the T = Drude-Boltzmann semi- 
classical conductivity of a 2D (or 3D) "metallic" system 
vary as a function of its carrier density in the presence of 
background screened Coulomb disorder being the main 
resistive scattering mechanism? 

In the rest of this article, we provide a detailed theo- 
retical answer to the question posed in the last sentence 
above. 



II. BACKGROUND 

The r = conductivity cr(n) of a 2D metallic system 
is typically written a o^'^° 



cr(n) 



'-D{Ef)t{Ef), 



(1) 



in the Boltzmann theory, where Ep^ vp are respectively 
the Fermi energy and the Fermi velocity, D{Ep) is the 
carrier density of states at the Fermi surface, and r is 
the so-called relaxation time (or the scattering time) for 



the relevant resistive scattering mechanism. (The factor 
2 in the denominator of Eq. ([T]) is replaced by 3 for 3D 
systems.) It is assumed that vp, Ep, and D are known 
as a function of the carrier density n from the relevant 
band structure information (and we will consider only 
the parabolic and the linear band approximation with 
the linear approximation used for obtaining our graphene 
transport results). The whole theory then boils down to 
a calculation of the transport relaxation time r at the 
Fermi surface using the appropriate microscopic scatter- 
ing mechanism (which we describe in details in section 
III of our paper) . 

For parabolic bands with i?(k) — h?k^/2m we can 
write Ep = mv^/2, where m is the carrier effective mass, 
and Eq. ([T]) simply reduces to the celebrated Drude- 
Boltzmann transport formula: 



and 



a{n) 



(2) 



where r = t{Ef) and n is the relevant (2D or 3D) 
carrier density. For graphene, Eq. ([2]) requires a slight 
modification^ which we will discuss later when we come 
to describing our graphene results. From this point on, 
unless otherwise stated, our explicit equations and for- 
mula are given for parabolic band 2D systems — we will 
consider the very special linear band case of graphene 
separately at the appropriate juncture. Most of our re- 
sults focus on 2D semiconductor systems, specifically 2D 
GaAs electron or hole quantum wells where the parabolic 
approximation applies well. We will point out the corre- 
sponding 3D parabolic band analytic results as appropri- 
ate, concentrating on presenting the equations and for- 
mulas mainly for 2D quantum well systems since most of 
our presented results are for 2D semiconductor systems, 
where the corresponding experiments are feasible. 

We note that one can, instead of discussing the con- 
ductivity (J, equivalently discuss the relaxation time r or 
the resistivity p = a^^ or the mobility /i, which is defined 
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Obtaining the density dependent conductivity a{n) now 
becomes equivalent to obtaining the density dependence 
of r or ^ since a ex nfi (or nr). We write: 



j3{n) = a{n) + 1 = 



dhic 
dhT 



n 



(7) 



We note that, in general, the mobility exponent a = 
a^ and the relaxation rate exponent ar may be unequal 
(e.g., graphene) except that for parabolic bands we have 
ctfj. — cur — ck. But the relation /S — afj, + l — a + 1 
always applies. The main goal of the current work is to 
calculate the exponent a{n) for various specified disorder 
mechanism, and discuss/contrast how the density scaling 
exponent depends on the type of disorder controlling the 
resistive scattering mechanism. 

We will obtain a(n) analytically in various limiting sit- 
uations in our theoretical study. But the main technique 
would be to obtain /i(n) or a{n) numerically from the 
Boltzmann theory and then to calculate a(n) or /3{n). 
Experimentally, the way to estimate the exponent a{n) 
is to plot the measured low-temperature /x(n) against n, 
and then to estimate a{n) by carrying out the logarith- 
mic differentiation. In general, a{n) will depend strongly 
on the background disorder, and will vary smoothly as a 
function of carrier density as different scattering mech- 
anisms are operational in different density regimes and 
as the background Coulomb disorder is screened differ- 
ently at different carrier density through the variation in 
qTp/2kp. An important surprising finding of our theory 
is an interesting nonmonotonic variation in the scaling 
exponent a{n) as a function of carrier density. 

It is important to point out at this stage that graphene, 
because of its linear dispersion with a constant Fermi ve- 
locity {vp = Wo), does not obey^ the exponent scaling re- 
lation (i = Ut -\-\ connecting the density scaling between 
conductivity (/?) and relaxation rate (a^). Instead, for 



graphene*, we have /3 



1/2 if a-r is defined by 



T (X n 



fj, (X n 



(4) 



T~ ~ n"^ which follows from the constancy of vp and 
the fact that D{Ep) ex Ep ex y/n in graphene. Putting 
D{Ep) (X kp (X y/n in Eq. ([T7]), we get for graphene 
a{n) ~ \/nT(n), i.e., 13 — ar + 1/2. Thus, whereas in 
ordinary parabolic systems the exponent a {= a-r = a^) 
is the same for both mobility (/x) and the relaxation rate 
(t~^), in graphene, by virtue of its linear band disper- 
sion, a = a^ = a^ — 1/2, but the relationship between 
the conductivity exponent (/3) and the mobility exponent 
(a) is still given hy P — ol-\-1 since by definition a = ne/i. 



and therefore 



,/3=a+l 



In general, we can define the conductivity exponent /3, 
or equivalently the mobility exponent a through the re- 
lation: 



a{n) 



rfln/i 
dlnn' 



(5) III. MODEL AND THEORY 

The central quantity to obtain theoretically in the 
semiclassical Boltzmann transport theory in the relax- 
ation time r or equivalently the relaxation rate t~^, 
which is given by the following expression for 2D sys- 
(.") tems, within the leading order Born approximation, for 



carrier scattering at T = by disorderiSiiS: 



1 _ 27r,-^ 
7 ^ "fi~2^ 



Nl''\z)dz 



cPk' 



vi'Xiz) 



(2^)2 
X (1 - cos6lkk')^[-E(k) - E{\i!)]. 



(8) 

Here N^'{z) is the 3D density of random impurities of 
the 7-th type (in general, there could be several different 
types of impurities present in the system: near and far, 
2D or 3D, long- or short-range) with z being the direc- 
tion perpendicular to the plane of 2D system (which lies 
in the x-y plane); Vq{z) is the electron-impurity interac- 
tion (in the 2D wave vector space defined by q) . k, k' 
are the incoming and outgoing 2D carrier wave vectors 
involved in the scattering process with a scattering angle 
^kk' between them and k — k' being the scattering wave 
vector; i?(k) = H'^k^/2m is the carrier energy. Note that 
our disorder model assumes a random distribution of im- 
purities in the 2D x-y plane although it is easy to include 
correlations in the 2D impurity distribution if experiment 
indicates the importance of such correlations. 

Once the scattering potential Vq{z) is defined, the 
problem of calculating the 2D conductivity becomes sim- 
ply a question of evaluating the 4-dimensional integral 
given in Eq. (JH). Note that we are restricting to the 
T = case (or to low temperatures), otherwise a ther- 
mal average would be required in defining t~^ necessitat- 
ing a 5-dimensional integration. We note that although 
Eq. ([8|) applies only to 2D systems, a very slight mod- 
ification gives us the corresponding expressions for 3D 
systems and graphene, which we do not show here. We 
note here that for graphene there is a well-known addi- 
tional form factor of (1 -I- cos^kk') inside the integral in 
Eq. ([5]) arising from chirality^. 

It is worthwhile to point out here that the theoretical 
idea of a meaningful universal density scaling behavior of 
conductivity applies as a matter of principle only when 
the resistive scattering is dominated by a particular dis- 
order mechanism. If there are many different types of 
disorder (i.e., several 7- types) contributing equivalently 
to the resistivity, then the net resistivity will be given 
by the Matthiessen's rule: T^^(n) — J2y'''y^i^)' ^'^d 

P{n) = J^-yPii'Ti), i-e-, <J-^{n) = J^.y'^^^i''^)^ where p^, 
a-y, T-y are respectively the resistivity, conductivity, scat- 
tering time for the 7-th type of disorder arising from 
Nl' {z) in Eq. ([5]). In such a situation, unless one par- 
ticular type of disorder (i.e., one specific 7) dominates 
scattering, the resulting density dependence of the total 
a (or /i) will manifest complex crossover behavior arising 
from the combination of all different scattering processes 
contributing with different strength, and there will not 
be only universal density dependent scaling behavior of 
u{n) or ^{n). To exemplify this important point, we 
consider a strict 2D electron gas being scattered by three 
different types of disorder (i.e., 7 = 1, 2, 3) given by 
remote random charged impurities at a distance d from 
the 2D system, background charged impurities at the 2D 
layer, and zero-range disorder in the layer, each with their 



respective conductivity exponent /3^ with 7 = 1, 2, 3 re- 
spectively. We can formally write: 



-1/' \ -1 I -1 I -1 
a [n) = 0-1 + (72 +0-3 
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(9) 



where A^ oc n^^^ is the strength of each independent 
physical scattering process (i.e., n^^i, 71^,2, "-i,3 denote re- 
spectively the remote and background charged impurity 
density and the short-range defect density). If we now 
define the net conductivity exponent /3(n) through the 
usual /3 ~ dlncr/dlnn definition using the total conduc- 
tivity (T(n), then obviously /3(n) will be a complex (and 
opaque) function not only of /3i, /32, and /^s, but also of 
the disorder strength rii^i, 71^,2, '^■1,3- Thus, the extrac- 
tion of a universal density exponent /3 (or a) makes sense 
only when one scattering mechanism dominates — gener- 
ically, the conductivity /mobility exponent j3/a (= /?— 1) 
is nonuniversal and depends in a complex manner both 
on the individual scattering mechanism and the relative 
strengths of different operational scattering mechanisms. 
We consider below theoretically several different disor- 
der potentials which may be operational in real 2D and 
3D systems, obtaining asymptotic analytical expressions 
for the exponent a and j3 in the process as applicable. 



A. Zero-range disorder 

Zero-range disorder, Vq = Vq (a constant), corresponds 
to pure (5-function real space scatterers distributed ran- 
domly spatially. Without any loss of generality, we can 
drop the z-dependence of the disorder (since the electron- 
impurity interaction is spatially localized), and assume 
the 2D electron system to be of zero thickness in the 
z-direction interacting with the random zero-range scat- 
terers situated in the same plane. For 3D systems we of 
course assume the scatterers to be randomly distributed 
three dimensionally. 

For the constant V"q model potential it is straightfor- 
ward to do the momentum integration in Eq. ([SJ to ob- 
tain the following results for 2D and 3D semiconductor 
systems and graphene: 



n° 2D, 

n^i/^ 3D, 



n 



-1/2 



graphene. 



(10a) 
(10b) 
(10c) 



We note that in graphene the mobility exponent {a = a^) 
differs from the relaxation rate exponent a-r by a^ — 
a-r — 1/2. In obtaining the density dependence of the 
relaxation time above we use the standard expressions for 
kp and Ep for parabolic 2D, 3D systems, and graphene: 

kp = (27rn)i/2; Ep = h^Tm/m : 2D, (11a) 

kp = (37r2n)i/3; Ep = h^{iTi'^mf/^ /2m : 3D, (lib) 
kp = (7rn)^/^; Ep = hvoiTrny^'^ : graphene. (lie) 



We use n throughout to denote the relevant 2D or 3D 
carrier density of the "metal", and Vq in Eq. ((Ht) is 
the constant graphene Fermi velocity defining its linear 
energy-momentum relationship, E — hv^k (wo ~ fO^ 
cm/s). We assume a spin degeneracy of 2 throughout 
and an additional valley degeneracy of 2 for graphene in 
defining kp and Ep. 

For zero-range (5-function disorder which is equivalent 
to assuming an uncorrelated white-noise disorder (often 
also called in the literature, "short-range disorder" , some- 
what misleadingly) , we therefore have (we define a = a^ 
always as the mobility exponent) 



(2D); -1/3 (3D); -1 (graphene), (12) 



and 



P^a+l^l (2D); 2/3 (3D); (graphene). (13) 

These results for 2D parabolic system and graphene are 
known^'W aj^j show that the conductivity grows linearly 
with carrier density in 2D systems and becomes a den- 
sity independent constant in graphene when transport is 
limited or dominated by zero-range white-noise type 5- 
function background disorder. The zero-range transport 
result for the 3D systems (metals or doped semiconduc- 
tors) does not appear to be as well-known (perhaps be- 
cause the density dependence of conductivity is not of 
much experimental interest in 3D systems as discussed 
in the Introduction of this paper) and shows surprisingly 
a sub-linear ^^ t?I^ increase in a{n) in 3D systems for 
(5-correlated short-range white-noise disorder. 

We now move on to the more interesting and experi- 
mentally more relevant Coulomb disorder, and discuss a 
number of disorder models pertaining to Coulomb disor- 
der in the next two subsections. 



B. Unscreened (long-range) Coulomb disorder 

In most semiconductor systems (2D or 3D), the dom- 
inant background disorder arises from quenched random 
charged impurities in the environment. Thus, the bare 
electron-impurity interaction is invariably the long-range 
1/r Coulomb interaction arising from the electric po- 
tential of the charged impurity. The charged impurity 
could be an intentional dopant impurity introduced to 
provide the doping necessary to create free carriers in 
the semiconductor or an unintentional charged impurity 
invariably present even in the cleanest semiconductor. In 
general, the Coulomb disorder from the random charged 
impurities should be screened by the free carriers them- 
selves so that the effective (screened) Coulomb disorder 
is short-ranged (this is the infra-red regularization nec- 
essary for handling the long range nature of Coulomb 
interaction). We consider the screened Coulomb disor- 
der in the next three sub-sections, focusing here on the 
unscreened bare Coulomb disorder for the sake of theo- 
retical completeness. 



The disorder potential Vq(£) is given by the following 
equation for the unscreened Coulomb interaction: 



^q(^) - 



nq 



.^q\A 



for 2D systems and graphene, and 



^q(^) 



47rZe2 



(14) 



(15) 



for 3D systems. Here q — |q| in Eqs. P^ and p^ is 
the appropriate 2D and 3D wave vector, Z is the atomic 
number of the charged impurity center with Ze being 
its charge (we consider Z = 1 throughout), and k, is the 
background static lattice dielectric constant. The coor- 
dinate z in Eq. ([Tl)) denotes the spatial separation of the 
charged impurity from the 2D confinement plane of the 
electron layer (taken to be located at z = in the x-y 
plane). We note that Eqs. p^ and p^ are simply the 
Fourier transform of the 1/r three-dimensional Coulomb 
interaction (with r = [x, y, z) a 3D space vector) in 2D 
and 3D systems. 

For 3D systems, unscreened Coulomb disorder leads to 
the following expression for the scattering rate: 



27rN, 



(fik' 

^2^ 



47re"^ 



Klk-k' 



X (l-cos0kk')'5(^fc-^fc'), (16) 

where all wave vectors are now three dimensional and Ni 
is the 3D impurity density. The momentum integration 
in Eq. p6)) is straightforward, and it is well-know n^ ^i^^ 
that the integral has a logarithmic divergence arising 
from the long-range nature of bare Coulomb interaction 
(i.e., an infra-red singularity). We get: 



(nln6) 



(17) 



with h 



^kl/qlp 



0, where kp — (37r^n)^/^, and 
qxF, which in principle is the 3D screening wave vector, 
goes to zero in the unscreened approximation, leading 
to a logarithmic divergence without the screening cut- 
off of the long-range Coulomb interaction. There are 
several ways of cutting off this long-range logarithmic 
Coulomb divergence (e.g., Conwell-Weiskoff approxima- 
tion or Brooks-Herring-Dingle approximation) which has 
been much discussed in the transport literature on doped 
semiconductorsiiii^,. Since our interest in mainly focused 
on 2D systems where the carrier density can be tuned 
continuously (in contrast to 3D systems), we do not fur- 
ther discuss the implications of Eq. pT)) for 3D doped 
semiconductors. 

For 2D systems, one must distinguish among (at least) 
three different kinds of Coulomb disorder: Random 2D 
charged impurities in the 2D plane of the carriers [i.e., 
z = in Eq. p4)) ]. random charged impurities in a 2D 
layer parallel to the 2D carriers with a separation d (i.e. 
z — d), and 3D random charged impurity centers in 
the background (where z varies over a region in space). 



We refer to these three situations as 2D near impuri- 
ties, remote impurities, and 3D impurities, respectively, 
throughout this paper. 

We first consider the 2D near impurity case with z = 
assuming unscreened impurity-electron Coulomb interac- 
tion in the 2D plane: 



V,iz = 0) 



2ne' 
Kq 



(18) 



This then gives (with n^ now as the 2D impurity density) : 



T-l = 



27rn,- 



(2^ 



27re^ 



Ik-k' 



X {l-co8e^^,)S{Ek-Ek' 



(19) 



The 2D integral in Eq. (TTOl . in contrast to the corre- 
sponding (divergent) 3D integral defined by Eq. P^ . is 
convergent: 

r-i - kp^ - n^^. (20) 

Thus a = 1 in 2D systems for the unscreened 2D im- 
purity case (while the corresponding 3D case is logarith- 
mically divergent). The exponent /? = a + 1 = 2 for 
the unscreened 2D Coulomb impurity in 2D "metallic" 
systems. 

Next we consider the 2D system with 2D remote im- 
purities (z = d 7^ 0). The relaxation rate [Eq. ([5])] is now 
given by 



27rn,- 



d^k' 

(2^ 



Klk-k'l 



xil-cos0kk,)SiEk-Ek'). (21) 

The integral in Eq. (|21l) can be rewritten in a dimension- 
less form 



kp 



-2xdo 



dx- 



VT 



~2 ' 



(22) 



where do = 2kpd is dimensionless. [We note that putting 
do = for d = gives us the 2D near impurity results 
of Eq. ((20| .] The asymptotic carrier density dependence 
{kp (X \/n) implied by Eq. ([22|) depends sensitively on 
whether kpd [= do in Eq. (|22|) ] is small [kpd <^ 1) or 
large {kpd^ 1). For kpd ^ 1, we get from Eq. 



kp 



and for kpd ^ 1, we get from Eq. 



kl 



-3/2 



(23) 



(24) 



Thus, a = 1 (3/2) for kpd < 1 (> 1) and therefore 
/3 = 2 (5/2) for kpd < 1 (» 1) for 2D carriers in the 
presence of remote Coulomb scatterers. 

Finally, we consider the 2D carrier system with a 3D 
random distribution of charged impurity centers. The 



integral for the 2D relaxation rate (Eq. ([8])) now becomes 
with qte -^ 



'?''"(S)^»""'"(g)' <^^' 



which has the same logarithmic divergence for the un- 
screened Coulomb disorder as the corresponding 3D case 
considered above in Eq. (|T7)) . but with a different 2D 
density exponent (^ n^''") from the corresponding 3D 
density exponent (~ n~^^^) in Eq. (fTT]) . Thus, 2D car- 
rier systems with unscreened 3D Coulomb disorder would 
have logarithmically divergent resistivity necessitating a 
length cut-off on the long range part of the bare Coulomb 
potential similar to the well-known situation for long- 
range bare Coulomb disorder in 3D systems (e.g. doped 
semiconductors) . 

It is important to emphasize our interesting theoretical 
finding that, although 3D unscreened Coulomb disorder 
leads to logarithmic long-range divergence in the resis- 
tivity of both 2D and 3D systems independent of the 
carrier density (i.e. t~^ diverges logarithmically), the 
corresponding situation for 2D carrier systems with 2D 
Coulomb disorder has no divergence and does not require 
any cut-off dependent infra-red regularization. It turns 
out that 2D metals with unscreened Coulomb disorder 
arising from random 2D charged impurities is perfectly 
well-defined within the Boltzmann transport theory. Of 
course, the unscreened Coulomb disorder model is not 
realistic and the calculated conductivity may not agree 
with the experimental data, but theoretically it is per- 
fectly well-defined with the density exponents a {/3) being 
1 (3/2) and 2 (5/2), respectively depending on whether 
kpd ^ 1 or 3> 1 where d is the location of the impurities 
with respect to the 2D carrier layer. 

One may wonder about the fundamental reason un- 
derlying the necessity for infra-red regularization for the 
3D unscreened case and not for the 2D unscreened case. 
This arises from the 3D bare Coulomb potential (~ l/f?^) 
being much more singular in the long wavelength q -^ 
limit than the corresponding 2D Coulomb potential (^ 
l/q). It turns out that this difference between 2D and 3D 
is sufhcient to make the 3D Coulomb disorder case infra- 
red divergent whereas the 2D case being non-divergent 
without infra-red regularization. 

We do not discuss here the unscreened Coulomb dis- 
order case for graphene since for graphene the density 
scaling of the conductivity is the same for both un- 
screened and screened Coulomb disorder (and we con- 
sider screened Coulomb disorder in the next subsection) 
by virtue of graphene screening wave vector qxp (oc kp) 
being proportional to the Fermi wave vector leading to 
both screened and unscreened Coulomb disorder having 
the same carrier density dependence in the conductivity. 



C. Screened Coulomb Disorder 

This is the most reahstic (as weh as reasonably compu- 
tationally tractable) model for calculating the resistivity 
due to Coulomb scattering within the Boltzmann trans- 
port theory. The basic idea is to use the appropriately 
screened Coulomb potential in Eq. ([8]) for calculating the 
relaxation rate. We use the conceptually simplest (static) 
random phase approximation (RPA) for carrier screening 
of the long-range electron-impurity Coulomb interaction. 
This means that the screening functions in 2D and 3D 
parabolic systems are the static dielectric functions first 
calculated by Stern~. and Linhard-*^^, respectively. For 
graphene, we use the dielectric screening function first 
calculated in Ref. [15|. 

The relevant screened Coulomb disorder potential is 
given by: 



Va 






(26) 



where Vq is the long-range Coulomb potential and e{q) is 
the appropriate static RPA dielectric function. In gen- 
eral, the screened Coulomb potential can be rewritten 
as 



2^e2 

K{q + QTEJ 



for 2D systems and graphene, and 



(27) 



Va 



iire^ 



K{q'^ + q^p) ' 



(28) 



for 3D systems. The screening wave vector qrF, some- 
times referred to as the Thomas-Fermi wave vector is 
given by the following expression (obtained in a straight- 
forward manner from the corresponding static polariz- 
ability function or the dielectric function in Refs. KLSMla ): 



Qtf 



2me 



2D, 



1/2 



qTF = f ^^~| 3D, 



qTF 



Ae^kp 
kHvo 



graphene. 



(29a) 
(29b) 
(29c) 



We have used a valley degeneracy of 1 (2) for 2D/3D 
(graphene) systems in Eqs. ([29]) . and chosen a spin de- 
generacy of 2. We note the (well-known) results that 
in the 2D parabolic electron system the screening wave 
vector is a constant, whereas in graphene (3D parabolic 
system) it is proportional to kp {y/kp). 

It may be worthwhile to discuss the dimensionless 
screening strength characterized by the parameter q^ = 



qxp/'ikp, which is given by 
2me^ 



1/2 



qs = 



Kh?kp 

4me^ 
TTuh^kp 
4e2 



2D, 



1/2 



-1/6 3D, 







kHvq 



graphene. 



(30a) 
(30b) 
(30c) 



We note the curious (albeit well-established) result 
that the dimensionless screening strength increases (very 
slowly in 3D ^ n^^^^) with decreasing density (~ n"^/^ 
in 2D) except in graphene where qxp on kp leading to a 
density- independent qs- Thus, the strong (weak) screen- 
ing limit with qg = qTp/2kp ^ 1 (<C 1) is reached at low 
(high) carrier density in 2D and 3D metallic systems. 
This peculiar density-dependence of screening has quali- 
tative repercussion for the transport exponents a and /3 
as a function of density. 

We now rewrite Eq. ([8]) for the 2D relaxation rate t~^ 
in terms of screened Coulomb disorder obtaining; 



7rft3fc2, 



/2^ 



l22iqs,do), 



with 



h2{qs,do) 



-2xdo^2^^ 



(x + qs)^Vl - a;2 



(31) 



(32) 



for 2D carriers and 2D impurities (with impurity density 
Ui per unit area) with qg = qj-p/^kp and do — 2kpd 



(with z 



with 



d defining the impurity locations) 

2 



Sirh^kl, 
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^23 (9s 



hsiqs) = 



dx 



{qs + Vl - x'^) 



2^2 ' 



(33) 



(34) 



for 2D carriers and 3D impurities (with impurity density 
Ni per unit volume). We note that putting q^ — (i.e. 
no screening) immediately produces the density scaling 
exponents obtained in the last subsection. 

For 3D carriers with (obviously) 3D random charged 
impurity distribution, we get 



where 



-^33(^5 



Niva 
STTh^kj, 



o'^a 



hsiqs), 



l-x 



^q's)' 



(35) 



(36) 



Again, putting q^ = in the 3D result above produces 
the logarithmically divergent relaxation rate discussed in 
Sec. IIIB for the 3D unscreened Coulomb disorder. 



Finally, for graphene we write down the relaxation 
rates for scattering by screened Coulomb disorder arising 
2D and 3D charged impurity distributions respectively by 
following the standard references^iii^: 



where 



nfi^vokF 



Ig2 ^ I dx 



27re" 



lG2iqs,do), 



'VT 



(x + qs)^ 



^ — 2xdo 



(37) 



(38) 



for 2D charged impurities located a distance d (with do 
2kpd) from the graphene plane, and for 3D disorder: 



N, 



nh'^vok'i, 



27re" 



iGsiQs), (39) 



with 



/g3 = / dx- 



{x + qs 



(40) 



and Ni being the 3D impurity density. In the 3D 
Coulomb disorder case for graphene, the random charged 
impurities are distributed in the graphene substrate with 
a uniform random 3D distribution with a 3D impurity 
density of Ni. 

For our numerical calculations of transport properties 
(to be presented in the next section) which would focus 
entirely on 2D systems with 2D impurities (both near and 
remote) , we will include the realistic width of the quan- 
tum well through a subband form-factor modifying the 
Coulomb matrix element arising from the finite thick- 
ness of the quantum well in the z-direction. This is a 
nonessential complication (making our numerical results 
compatible with and comparable with the experimental 
low-temperature transport data in GaAs quantum wells) 
which does not affect our theoretical conclusions about 
the carrier density scaling of the 2D transport properties 
since the quasi-2D quantum well form factor is indepen- 
dent of the carrier density in the leading order. 

Below we obtain the asymptotic density exponents 
(based on the results given above) a and /3 for screened 
Coulomb disorder in the strong (q^ 3> 1) and weak 
(qs <S^ 1) screening situations considering both near and 
remote 2D impurities and 3D impurities. 



D. Strong-screening (qs 3> 1) and weak-screening 
{qs < 1) limits 

It is straightforward to carry out the asymptotic 
expansions of the various integrals in Eqs. ([5Tt - pn|) 
to obtain the strong-screening (qs ^ 1) and the weak- 
screening {qs <C 1) limiting behaviors of the relaxation 
rate t~^ for the different systems under consideration. 
Remembering that r^^ ~ n~°' and (3 = a + 1 {a + 1/2 



for graphene) we get the following results by taking 
(7s ^ 1 and qs -€1 I limits of Eqs. ((3T |) -(|40 | in various 
situations. 

(i) 2D carriers with 2D impurities: 

For 2fci?(i <C 1 (i.e. near impurities) 

a — {P = 1) for strong screening {qTF ^ 2kp) 

a = 1 {f3 = 2) for weak screening {qxr ^ 2kp). 

For 2kpd ^ 1 (i.e. remote impurities) 

a = 3/2 (/3 = 5/2) for both weak {qrp < 2kp) and 

strong {qxp ^ 2kp) screening. 

(ii) 2D carriers with 3D impurities: 

a — 1/2 (/3 = 3/2) for strong screening {qpp ^ 2kp) 

a = 3/2 (/3 = 5/2) for weak screening {qpp <C 2kp). 

(iii) 3D carriers with 3D impurities: 

a — 1/3 (/3 — 4/3) for strong screening {qrp » 2kp) 

a — 1 {13 — 2) for weak screening {qxp ^ 2kp). 

(iv) Graphene with 2D impurities (remembering 
a = Ot - 1/2 = a^): 

For 2kpd <C 1 (i.e. near impurities) a,- — 1/2; a = 
(/3 = 1) for both strong and weak screening. 
For 2kpd ^ 1 (i.e., remote impurities) ar — 1] a — 1/2 
(/3 — 3/2) for both strong and weak screening. 
We note that for graphene a = afj. = ar — 1/2 (and 
/3 = a -I- 1) by virtue of its linear dispersion. Also, 
for graphene strongly screened, weakly screened and 
unscreened Coulomb disorder manifest the same density 
exponent in the conductivity and the mobility since 
qxp oc kp, and thus qs is density independent. (Of 
course, the actual numerical values of the conductiv- 
ity and the mobility are very different in the three 
approximations except for having the same power law 
dependence on the carrier density depending only on 
whether the random 2D charged impurities are near or 
far.) 

(v) Graphene with screened 3D impurities: q;t- = 1; 
a — 1/2 (/? = 3/2) for both strong and weak screen- 
ing. 

We note that graphene with unscreened 3D Coulomb dis- 
order {qs = 0) has the same logarithmic divergence in 
the relaxation rate t~^ (and hence in the resistivity) for 
all carrier densities as in the corresponding 2D and 3D 
parabolic electron systems for unscreened 3D disorder. 
The unscreened 3D Coulomb disorder is thus unphysical, 
always necessitating an infra-red regularization as was al- 
ready realized in the 1950s in the context of 3D doped 
semiconductor transpor t ^^i^^ . 

In Table I we summarize our asymptotic analytic find- 
ings for the density scaling of conductivity and mobility 
for 2D, 3D parabolic systems and 2D graphene in vari- 
ous limiting situations involving Coulomb disorder aris- 
ing from random charged impurities in the environment. 

In Table I above we provide the asymptotic density 
scaling exponent (a) for the carrier mobility (/z ~ n") in 
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2D 


3D 


Graphene 


2D Coulomb disorder with near 
impurities {2kFd ^ 1) 


strong screening 
[qs » 1) 





N/A 





weak screening 

{is « 1) 


1 


N/A 





unscreened (g^ = 
0) 


1 


N/A 





2D Coulomb disorder with 
remote impurities {2kpd ^ 1) 


strong screening 
(9s » 1) 


3/2 


N/A 


1/2 


weak screening 
[qs « 1) 


3/2 


N/A 


1/2 


unscreened (g^ = 
0) 


3/2 


N/A 


1/2 


3D Coulomb disorder with 3D 
impurity distribution 


strong screening 


1/2 


1/3 


1/2 


weak screening 


3/2 


1 


1/2 


unscreened 


log-divergent 


log-divergent 


log-divergent 


zero-range disorder with (5-function 
impurities 


concept of 
screening inap- 
plicable here 





-1/3 


-1 



Table I. Asymptotic values of a. By considering the Coulomb disorder arising from random charged impurities in the envi- 
ronment the asymptotic density scaling exponent [a) for the carrier mobility (/i ~ n") in various limiting situations and for 
various types of disorder is given. The corresponding conductivity exponent (/3) with a ~ n^ is given \yy P — a-\- \. 



various limiting situations and for various types of disor- 
der. The corresponding conductivity exponent (/?) with 
a ^ nP is given by /3 = a -f 1. 



{Ae^kp I nhva) >• 2fci?: graphene, i.e. 



IV. RESULTS 



In this section we present our numerical results for the 
T = density dependent transport properties of 2D sys- 
tems within the Boltzmann transport theory in order to 
obtain the full density dependence of the exponent a{n) 
for mobility and equivalently the exponent /3(n) for con- 
ductivity, obtaining in the process the density regimes 
where the analytical asymptotic exponents obtaining in 
the last section apply. The reason we focus on the 2D car- 
rier system is that it is the most convenient system for 
the experimental investigation of the density dependence 
of transport properties. In 3D semiconductor systems, 
the carrier density cannot be continuously tuned as it 
can be in 2D systems. 

We first note that the strong (weak) screening con- 
dition implies low (high) values of carrier density in 
the system. Using Eq. pT|) and Eq. (f^ . we find: 
q^ = qrpp/2kF > 1 implies 2me^ / kH^ > 2(27rn)^/2. 
2D, (4meV^K;i2)i/2(3^2^)i/6 ^ 2(3^2^)1/3; 3D, 



( 



2me 
in \ Kh? 

2e^ 



>n 2D, 



>n 3D, 



^ va graphene. 



(41a) 
(41b) 
(41c) 



Eq. (|4T|) define the low-density regime where the 
strong screening condition would be satisfied except for 
graphene which has qs independent of carrier density 
since qrp oc kp. From Eq. (I41[) we conclude that strong 
(weak) screening situation (within RPA) for T = trans- 
port properties would be satisfied under the following 
conditions for the different systems under consideration: 

2 



n< (>) 



n< (>) 



X 1.14 X Wcm,-^ 2D, (42a) 



X 7.4 X 10 



21, 



3D, 



(42b) 
(42c) 



K <C (^) 4.4 graphene. 

Using the band effective mass for GaAs electrons [m = 
O.OTrrie) and holes (m = 0.4me), we get (using k = 13 
for GaAs-AlGaAs quantum wells) 

n < (>) 3.3 X 10"cm"2 for 2D n - GaAs, (43a) 

n < (>) 9.5 X 10"cm"2 for 2D p - GaAs, (43b) 

and 

n < (>) 1.5 X lO^^cm^^ for 3D n - GaAs, (44a) 

n < (>) 2.8 X IQi^cm-^ for 3D p - GaAs. (44b) 
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Figure 1. (color online) Calculated mobility, iMi, as a func- 
tion of density n of n-GaAs quantum well with a well width 
a = 200 A for various d =0, 200, 400, 600, 800 A(from bot- 
tom to top). Here rii is the 2D random charged impurity 
density located a distance d away from the quantum well. 
The mobility is calculated at T = K. 
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Figure 2. (color online) Calculated mobility exponent a in 
fj, oc n"' (i.e., a — dlnfi/dlnn) for the corresponding /i(n) in 

Fig.ra 



We note (again) that in graphene (as Table I indicates) 
the exponents a, /3 do not depend on weak/strong screen- 
ing or on the carrier density. The density range for GaAs- 
based 2D systems, which we would consider numerically 
in this section, is the experimentally relevant 10^ — 10^^ 
cm^'^ density range in high-mobility GaAs-Ala;Gai_a;As 
2D quantum well structures, and thus for 2D electron 
systems (2D n-GaAs) the crossover from the low-density 
strong-screening to high-density weak-screening behav- 
ior may be experimentally relevant. On the other hand, 
for 2D p-GaAs hole quantum wells, the crossover density 
(~ 10"'^* cm~^) is too high to be relevant experimentally, 
and thus the laboratory 2D hole systems are likely to be 
always in the strongly screened situation. 

The 2D transport exponents a and /3 depend on an 
additional dimensionless parameter (i.e. in addition to 
Qs = qTE/'^kp) do = 2kpd, which depends both on the 
carrier density through kp ~ \/n and on the separation 
{d) of the random charged impurities from the 2D sys- 
tem. This dependence on kpd, the dimensionless sepa- 
ration of the impurities from the carriers in the z direc- 
tion (which has no analog in the corresponding 3D dis- 
ordered case), is experimentally always relevant in high- 
mobility 2D semiconductor structures for three reasons: 
(1) in high- mobility modulation-doped 2D quantum well 
structures, scattering by remote dopants (which are in- 
troduced intentionally at a known distance d from the 
quantum well) is always present; (2) in principle, it is 
always possible to place ionized impurities at a set dis- 
tance d from the 2D quantum well using the computer- 
controlled MBE growth technique which is used to pro- 



duce high-quality semiconductor quantum wells to begin 
with; (3) in realistic experimental samples with a finite 
quasi-2D thickness, d is always finite. 



In general, the experimental 2D mobility/conductivity 
could be limited by various types of Coulomb disorder—: 
near/far 2D Coulomb disorder {2kpd <^ 1/ ^ 1): 3D 
Coulomb disorder in the background. According to Table 
I, the asymptotic density dependence in each case should 
be different. In Figs. [T] and [2] we show our directly numer- 
ically calculated n-2D GaAs mobility for quantum well 
electrons assuming 2D random charged impurity scat- 
tering from quenched point scattering centers located a 
distance d away from the quantum well (including the 
d — case). The numerically calculated mobility expo- 
nent a{n) = dlnfi/dlnn is shown in Fig. [2] for the cor- 
responding /i(n) shown in Fig. [T]- we note that n (x kp, 
and therefore the parameter kpd ^ \fn for a fixed value 
of d. The results shown in Figs.[T]and[5]correspond to the 
zero-temperature case, but are indistinguishable from the 
corresponding low-temperature results for T — i^^niK 
(which we have verified explicitly). We note that the nu- 
merical results presented in Figs. [1] and [5] are realistic 
(as are all other numerical results shown in this paper) 
in the sense that they are obtained from the full numer- 
ical integration of the Boltzmann theory expression for 
the relaxation rate [as given in Eq. (jS]) in section III] 
with the additional sophistication of including the finite 
thickness of the quantum well through the finite well- 
thickness form-factors /i((7) and /(g) which modify the 
Coulomb disorder matrix element (i.e. \y\^ — )■ |l^P/i(g)) 
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Figure 3. (color online) Calculated mobility of n-GaAs quan- 
tum well (a = 200A) with unintentional charged impurities. 
Here the 2D impurities {ni2D = 10^ cm~^) are at the GaAs- 
AlGaAs interface whereas there are two different types of 
3D Coulomb disorder: inside the GaAs well [rii^n — 10^'' 
cm^^ (GaAs)] and inside the barrier AlGaAs [tusd = 10^^ 
cm-='(AlGaAs)]. 



and qTF — >■ iTpfiq), respectively, and they are given by 
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where a refers to the quantum well width taken to be 
200 A for the results in Figs. [Hand El We note that the 
form- factor fi{q) simply reduces the Coulomb impurity 
potential from its q^^ behavior by a g-dependent (but 
density-independent) function determined by the thick- 
ness a of the GaAs quantum well. The form- factor f{q) 
reduces the 2D screening through the modification of the 
electron-electron interaction due to the finite thickness of 
the quantum well. Since the quantum well form factors 
do not depend on the carrier density in the leading or- 
der, this quasi-2D form-factor effect does not in any way 
modify the asymptotic exponents a and (3 given in Table 
I (and theoretically defined in the lase section), but the 
form factors do modify the actual calculated values of the 
mobility/conductivity/resistivity, making them more re- 
alistic, and therefore any comparison with experimental 
transport data necessitates the inclusion of the quasi-2D 
form factor effect. All our numerical results for 2D elec- 
trons and holes in GaAs quantum wells presented in this 
paper include the realistic quantum well form factors in 
the theory taking into account the finite well width ef- 
fect both in the electron-impurity interaction and in the 
electron screening [with qTF being modified to qxpfiq)]- 
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Figure 4. (color online) Calculated mobility exponent a in 
fj, (X n" (i.e., a — dlnfi/dlnn) for the corresponding /i(n) in 
Fig. [3] 



In Fig. [T] we show (we actually show the calculated fxrii 
since /i ^ l/rii) our calculated mobility /i for 2D n-GaAs 
as a function of carrier density n for different impurity 
locations (d) whereas in Fig. [5] we show the corresponding 
mobility exponent a{n) — dln/i/dlnn. The asymptotic 
high-density results, a — >■ 3/2 for n -^ oo (i.e., qxp <C 
2kp) and kpd^ 1 as given in Table I, is clearly obeyed in 
all cases with a « 3/2 for larger d and n values satisfying 
2kFd > 1. For 2D GaAs systems: 



2kpd. 



5dVn, 



(46) 



where d is measured in 1000 A units and n in units of 
10^° cm-2. Thus even for d = lOOA, 2kFd > 1 already 
for n > 2 X 10^" cm^^, and thus the 2/ci?d :$> 1 condition 
is quickly reached at higher carrier density in all 2D GaAs 
systems for any type of relevant background 2D Coulomb 
disorder (since typically, n « 10® cm~^ is a lower limit 
for the achievable carrier density in 2D semiconductor 
systems). One may wonder if the d = situation corre- 
sponds to the 2kFd = situation considered in Table I. 
This is certainly true for the strict 2D limit (i.e. a — 
limit of the quantum well). But for any finite value of a, 
the d = impurity location only refers to the distance 
of the 2D charged impurities from the GaAs- AlGaAs in- 
terface, and thus the average impurity separation from 
the electrons is always finite except in the a — >■ limit. 
Thus even the d — case in our theoretical calculation 
has an effective finite value of do — 2k pd because of the 
finite layer thickness effect (e.g., d sa a/2 effectively in 
the d —^ limit). 

The most important conclusions from Figs. [T] and [2] are: 
(i) The 2kpd ^ 1 condition dominates the mobility ex- 
ponent except for rather low mobility samples with very 
small values of d; (ii) even for the nominal d = case 
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in Fig. 1 (where in the strict 2D case, a ^ 1 always) a 
eventually approaches the asymptotic a — >■ 3/2 value for 
n ^ 10^^ cm~^ because of the finite layer thickness effect 
(i.e., finite a); (iv) for scattering purely by very remote 
dopants {d > 500 A), the mobility exponent a > 1 al- 
ways because 2kFd ^ 1 condition is always satisfied; (v) 
the low density limit {2k pd ^ 1, qs ^ 1), where a — >■ 
according to Table I, would be achieved in 2D n-GaAs 
systems only for n <C 10^ cm~^, and in all realistic situ- 
ations, a > 0.5 always as long as transport is dominated 
by 2D Coulomb disorder. 

In Fig. 13] and m we show (again for a = 200 A) our cal- 
culated mobility in the presence of both 2D and 3D (un- 
intentional background) disorder neglecting remote scat- 
tering effects (assuming the intentional remote dopants 
to be too far, d > 1000 A, for them to have any quanti- 
tative effects, as would apply to gated undoped HIGFET 
structures or to extreme high-mobility modulation doped 
structures where the intentional dopants are placed very 
far away). In Figs. [3] and IH the 2D impurities (71.^20) are 
put right at the GaAs-AlGaAs interface whereas there 
are two different types of 3D Coulomb disorder: inside 
the GaAs well {riisD (GaAs)] and inside the barrier Al- 
GaAs [ni3D (AlGaAs)]. Again the corresponding criti- 
cal exponents increase with carrier density, approaching 
a = 3/2 for high density consistent with the analytic 
theory. When the dominant background disorder is that 
arising from riisD (GaAs), i.e., unintentional background 
impurities in the well itself, typically a w 0.5 — 0.8, which 
is in between weak and strong screening situation. 

In Figs. [5] and [SI we show our 2D p-GaAs results for 
hole-doped high-mobility GaAs quantum wells. Every- 
thing in Figs. [5] and |6] is identical except for using two 
different hole effective mass values: nih — 0.3mo (Fig. [5]) 
and 0.4too (FiglS]) since the hole effective mass in GaAs 
quantum wells is somewhat uncertain". The precise 
value of rrih affects qtf oc to, and thus determines the 
value of qs = qxp/'^kp, leading to some difference be- 
tween the results in Figs. [S] and H) For the holes, we 
show both the individual mobility and exponent for each 
scattering process as well as the total mobility and to- 
tal exponent obtained by adding the two resistivities (or 
equivalently, the two scattering rates arising from the two 
scattering mechanisms). We deliberately refrain from 
showing the hole results (Figs. [S]and[n]) in the same for- 
mat as the electron results (Figs. [1] and ^ since they 
would all look identical except for some changes in the 
numbers. 

In Figs.[S)/|n](a) and (c) we show the calculated 2D hole 
mobility /i(n) limited by 3D Coulomb scattering (rii^o) 
and 2D remote Coulomb scattering (n^) for a = 200 A 
with the only difference being rii = 8(4) x 10^^ cm~^ 
respectively in Fig. ^^ (a) (c), showing explicitly the 
quantitative importance of 2D remote Coulomb scatter- 
ing vis a vis 3D Coulomb scattering. The corresponding 
critical exponent a(n), shown in Fig. [S|ini (b) and (d) re- 
spectively, is completely consistent with Table I with a 
for remote scattering quickly reaching the asymptotic un- 
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Figure 5. (color online) (a) Hole mobility as a function of 
hole density (n) of p-GaAs quantum well with a width a — 
200 A and the hole mass of m = 0.3. Here the background 
unintentional charged impurities with a density Uiso = 3 x 
lO^^cm""^ are located inside the quantum well and 2D remote 
charged impurities with a density n^ = 8 x lO^^cm"'^ are 
located at d = 150 A from the interface. In this figure the 
black solid curve indicates the total mobility and blue dot- 
dashed (red dashed) curve indicates the mobility limited by 
only background scattering (remote charged scattering), (b) 
The exponents (a) for the corresponding mobilities in (a), (c) 
Hole mobility with the same parameters of Fig. [5ja) except 
the remote charged impurity density n^ = 4 x lO^^cm"^. (d) 
The exponents (a) of corresponding mobilities of (c). 



screened value of 3/2 as 2kFd ^ 1 and a for 3D Coulomb 
disorder increasing slowly from the very strongly screened 

to 
where 



low density situation {asn ^ 0.1 for n < 2 x 10^ cm 

10^2 en 



a ~ 1 for very high hole density (n 
screening weakens. 

The total exponent in Figs. [5] and [6] shows complicated 
non-monotonicity as a function of carrier density since 
the low-density (high-density) situation is more strongly 
affected by remote 2D (background 3D) Coulomb scatter- 
ing and the density dependent crossover between the two 
scattering regimes is completely nonuniversal depending 
precisely on the relative amounts of 2D and 3D Coulomb 
disorder (i.e., on ni^D^ Ui, and d). The only concrete 
statement one can make is that a{n) increases at low car- 
rier density generally to a value larger than unity whereas 
at intermediate to high density it tends to stay below 
unity. Again, all of these results are completely consis- 
tent with the asymptotic exponents given in Table I (as 
long as various scattering mechanisms with different ex- 
ponents are combined together). 

So far we have discussed our numerical results of 
Figs. [T]— [B] in terms of their consistency with the the- 
oretically analytically obtained critical exponents given 
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Figure 6. (color online) (a) Hole mobility as a function of 
hole density (n) of p-GaAs quantum well with a width a = 
200 A and the hole mass of m = 0.4. Here the background 
unintentional charged impurities with a density nisu = 3 x 
lO'^^ cm~^ are located inside the quantum well and 2D remote 
charged impurities with a density n^ = 8 x 10^^ cm"^ are 
located at d = 800 A from the interface. In this figure the 
black solid curve indicates the total mobility and blue dot- 
dashed (red dashed) curve indicates the mobility limited by 
only background scattering (remote charged scattering), (b) 
The exponents (q) for the corresponding mobilities in (a), (c) 
Hole mobility with the same parameters of Fig. [5]^ a) except 
the remote charged impurity density n^ = 4 x lO^^cm"^. (d) 
The exponents (a) of corresponding mobilities of (c). 



in Table I with the mobihty exponent (/u ~ n"'"') a 
showing the expected behavior in the asymptotic density 
regimes of 2kpd ^ 1 (<C 1) and Qs ^ 1 {<^ 1) as the case 
may be. Remote scattering by 2D ionized dopants dom- 
inates transport at low density {2kpd < 1) crossing over 
to background impurity scattering dominated regime at 
higher density, leading to a > 1 (< 1) at low (high) 
density. 

Now we discuss perhaps the most interesting aspect of 
our numerical results in Figs. [T]— El which appears to be 
at odds with our asymptotic theoretical analysis of sec- 
tion III (and table I). This is the intriguing result that 
a{n) due to remote 2D impurity scattering can actually 
exceed the asymptotic unscreened 2D Coulomb scatter- 
ing value of a = 3/2. It is clear in Figs. [H H E^b) , EJd) , 
inib), [BJd) that there is a shallow maximum in a{n) at 
some intermediate density where the numerically calcu- 
lated mobility exponent a > 3/2 — 1.5 (and therefore the 
corresponding conductivity exponent /3 > 2.5), which is 
remarkable since the unscreened Coulomb scattering (ap- 
plicable in the high-density regime defined by ikpd 3> 1 
and/or q^ >> 1) by remote impurities produces a = 3/2. 
This non-monotonic behavior of the exponent a{n) in 



the intermediate density regime (neither high- nor low- 
density asymptotic regime considered in Table I) with an 
exponent value larger than the corresponding unscreened 
Coulomb exponent 3/2 is unexpected and highly intrigu- 
ing. We provide a theoretical explanation for this intrigu- 
ing nonmonotonic behavior of remote Coulomb scatter- 
ing at intermediately density with a mobility exponent 
exceeding the unscreened value of 1.5 in the next sec- 
tion. 

Before concluding this section on numerical results, we 
provide, for the sake of completeness, our numerically 
calculated mobility exponent for 2D graphene transport 
under the remote Coulomb scattering situation. In Fig. [7] 
we show our calculated graphene mobility exponent a{n) 
obtained from the numerically calculated graphene mo- 
bility, a{n) — dhifi/dhin^ for various locations (d) of 
the 2D impurity layer in relation to the 2D graphene 
layer. We note that graphene is a strictly 2D system 
and hence there is no quasi-2D form factor correction. It 
is clear that a{n) goes asymptotically to the unscreened 
Coulomb value of 3/2 as n (and therefore kpd increases) 
except for the trivial d = case where a{n) — (i.e., 
(3 — 1) for all density as is already well-known in the 
literatur o^i^^ and is well- verified experimentally^^. 

The calculated graphene mobility exponent a{n) 
shown in Fig. [7] agrees completely with the analytical ex- 
ponent values given in Table I. (We mention again that 
in graphene the mobility exponent a^ = a and the relax- 
ation rate exponent ar differ with a = a^ = ar — 1/2 and 
/3 = a+l = afj_ + lhy definition whereas in 2D and 3D 
parabolic system, where fi oc t, a^ — ar ~ a = (3 — 1.) 
We note that in graphene, for impurities away from the 
2D graphene plane (i.e., d ^ 0), the asymptotic conduc- 
tivity exponent (3 for high carrier densities (kpd 3> 1) 
is 3/2 and thus a ex n^^'^ in graphene layers dominated 
by far away Coulomb impurities. An experimental veri- 
fication of such a cr ~ n^''^ behavior in graphene due to 
Coulomb scattering by remote impurities will be a direct 
verification of our theory. 



V. NONMONOTONICITY OF TRANSPORT 
SCALING 

We now theoretically discuss (and explain analytically) 
our surprising numerical finding in section IV, not antic- 
ipated at all in the asymptotic theory of section II or in 
any of the substantial earlier literature on 2D transport, 
that the density scaling exponent a {(3) of 2D mobil- 
ity (conductivity) has an intriguing nonmonotonicity as 
a function of carrier density in the intermediate density 
regime (in-between the asymptotic low and high density 
regimes discussed in section III and tabulated in Table 

I)- 

We first note that the nonmonotonicity in a{n) arises 
from the subtle fact that although the exponent a (or /3) 
depends only on one explicit external variable (namely, 
the carrier density n), it depends theoretically on two 
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Figure 7. (color online) (a) shows the calculated graphene 
mobility exponent otin) obtained from the numerically cal- 
culated graphene mobility, a(n) = dln/i/dlnn, for various 
location (d) of the 2D impurity layer in relation to the 2D 
graphene layer. In (b) the exponent for d = 30 nm shows a 
shallow local maximum at n ~ 10^^ cm~^. 



precisely a = 1 in the high density (g^ ^ 1) for d = 
as one would expect theoretically (we have verified this 
strict 2D limit with d = result explicitly numerically). 

For d 7^ (i.e., Ikpd ^ 0), however, the behavior of 
transport properties depends nontrivially on the variable 
do = 2kpd, and there is no obvious theoretical reason 
why the low-density (a = or 1 depending on strong or 
weak screening) and the high-density (a — 3/2 always) 
asymptotic limits must be the lower and upper bounds 
on the exponent. In fact, as our numerics show, and we 
establish theoretically below, a{n) does have a peak (ex- 
ceeding the unscreened a = 3/2 value) at an intermediate 
density around 2kFd^ 1. 

As shown in our numerical results presented in the last 
section, a{n) defined by /i ~ n"^"', can exceed the un- 
screened exponent value a = 3/2 reached in the asymp- 
totic high-density {2kpd :^ 1) limit. To verify whether 
this finding is a numerical artifact or real, it is sufficient 
to consider the unscreened remote 2D Coulomb disorder 
in the strict 2D electron layer limit where the free carriers 
and the charged impurities, both confined in infinite zero- 
thickness 2D layers in the x-y plane, are separated by a 
distance d in the z-direction. (This is the model explic- 
itly used in section III.) As n — > 0, a{n) — >■ 1 (unscreened 
2kpd <^ 1 limit) whereas as n -^ oo, a{n) — )• 3/2 (un- 
screened 2kpd 3> 1 limit). The zero-density limit would 
be modified to a = if screening is included in the theory, 
but we are interested here in the intermediate-density be- 
havior, not the zero-density regime. 

The scattering rate t~^ in this case is given by [see 
Eqs. dSJ) and (|5^]: 



C 



1 g — b2;Vn 



VI 



dx. 



(47) 



where we have shown the explicit density (n)-dependence 
everywhere (C and b are unimportant carrier density in- 
dependent constants for our purpose) and have used the 
fact that kp ^ y^. We rewrite Eq. (|47p as. 



C 



/(n). 



(48) 



where 



independent dimensionless variables g^ = qpp I2kp and 
do = 2kpd since in reality there are two independent ex- 
ternal variables in the problem: carrier density (n) and 
the impurity location (d). The dependence on two in- 
dependent variables is the key feature allowing for the 
presence of nonmonotonicity in a(n) as well as its max- 
imum possible value being larger than the unscreened 
exponent value a — ^ 3/2. Indeed in the strict 2D limit 
with d = (see, e.g., the graphene result in Fig. [T]), 
there is no maximum allowed in a(n). This is true for 
both graphene and 2D parabolic system - in graphene, 
a{n) = for all values of n in the d = limit whereas 
in strictly 2D parabolic system a(n) monotonically in- 
creases from a = in the low-density limit (g^ ^ 1) to 



I{n) 



-hx^/n 



vr^ 



^dx. 



(49) 



It is obvious that t~ is a monotonic function of n de- 
creasing continuously with increasing density with no 
extrema whatsoever since both 1/n and exp(— foxy^) 
decrease with increasing n continuously. This can be 
easily checked explicitly by showing that the equation 
dT~^{n)/dn — has no solution. The monotonic de- 
crease of r~^(n) with n simply implies that /x(n) oc T(n) 
increases monotonically with increasing density as is 
obvious from our numerical results in section IV: For 
Coulomb disorder, 2D mobility and conductivity always 
increase with increasing density monotonically. 
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To figure out whether a{n) = dlnfi/dlnn has non- 
monotonicity (or extrema) as a function of n, we must 
use /Lt(ri) oc r(n), and write 



d hi /i dlnfi 

a = — = n- 



d\nn 



dn 



n dl 
I dn' 



(50) 



where I{n) is the integral defined by Eq. (l49t . It is 
straightforward, but messy, to show that the condition 
dfi/dn = with a{n) defined by Eq. ([50]) has a solution 
at an intermediate value of n, and thus a{n) has an ex- 
tremum - it is still messier to show that d^a/d-n? < at 
this extremum so that a{n) has a maximum at an inter- 
mediate density as is found numerically in section IV. For 
our purpose, however, it is much easier to simply estab- 
lish that the function a{n) defined by Eq. (|50|) approaches 
the high density n ^ oo limit of a(n — >■ (X)) = 3/2 from 
above, thus definitively proving that the mobility expo- 
nent a exceeds 3/2 at some intermediate density (and 
thus must have a maximum in the < n < cxd or equiv- 
alently in the 1 ^ kpd^ 1 regime). 
As n — > 0, we have from Eq. (P^)) : 



by/n, 



leading to 



a{n ^ 0) = 1 + b^/n/TT. 
As n — > oo, we have: 



I{n) = 



1 



b^ 



1 



1 



leading to 



a{n ^ oo) = 3/2 + l/2b'^n. 



(51) 



(52) 



(53) 



(54) 



From Eqs. dni), (El, (SZl), and (gg) we have: b = 2V2^d, 
and thus 5 ex d is positive definite (except for the d = Q 
explicitly left out here). We, therefore, immediately 
conclude that a{n) approaches its asymptotic value of 
a — 3/2 for n — >■ oo from above with a{n — >■ oo) w 
I + Yg^^, and the leading possible correction to the 
exponent in the n — >■ oo limit is of 0{l/n) with a coef- 



ficient l/lQird 



ex 



We also find that the correction 



to the n — > value of a{n = 0) = 1 is of 0{\/n) with 
a positive coefficient of 2v27r(i/7r oc d. Thus, we now 
have a complete theoretical understanding of the intrigu- 
ing (and hitherto unexpected in the literature) numeri- 
cal finding in Sec. IV that, although the mobility ^{n) 
itself is a monotonically increasing function of increas- 
ing carrier density, its power law exponent a{n) shows 
a maximum (around 2fci?d ^ 1 in fact) approaching the 
asymptotic high-density [n — >■ oo; 2kFd S> 1) value of 
a = 3/2 from above, allowing a{n) to be larger than 1.5 
at some d-dependent value a > 1.5. 

One feature of our theory presented in Eqs. (|T7)) 
- (j54p is worth mentioning and comparing with the 
numerical results of Sec. IV. This is our finding in 



Eq. ((54| that the asymptotic high-density (n — >• oo) ex- 
ponent a{n — > oo) — 3/2 is approached from above as 
a{n -^ oo) = 3/2 + l/(167r(i^n), implying that the maxi- 
mum possible values of a, amax, scales approximately as 
(d^rimax)"^ ot {kpmdY, whcrc rimax and kpm are respec- 
tively the carrier density and the corresponding Fermi 
wave vector at the maximum. This implies that the maxi- 
mum value amax ~ 1-7 (~ 3/2-1- l/167r) is approximately 
independent of the value d and of the carrier effective 
mass with the value of the carrier density nmax (where 
the maximum occurs) scaling roughly as rimax ^ d^^. 
This strong prediction is approximately consistent with 
our numerical results - in fact, Qmax ~ 1-7 is clearly in- 
dependent of whether the system is a 2D electron or hole 
system and of the precise value of d. In fact, amax ^1-7 
being approximately independent of electrons/holes and 
the value of the separation distance d is a striking theo- 
retical result which is consistent with the full numerical 
results of Sec. IV. 

Because of the striking nature of our finding that 
Climax ~ 1-7 always (for remote impurity scattering) for 
2D electron/hole carrier systems, we carried out addi- 
tional numerical calculations using the realistic Boltz- 
mann theory (including both quasi-2D finite thickness 
and screening effects) for 2D n-GaAs wells of thickness 
a = 300 A (different from the case of a = 200 A used in 
Sec. IV) and incorporating both remote impurity scatter- 
ing with 2D impurity densities n^ and separation d and 
also (different values of Ud is used) with near impurity 
scattering with 2D impurity density Ui with d — Q (i.e. 
interface impurities). The calculated exponents for the 
individual scattering mechanisms a^ and a^ (for Ud and 
rii^ respectively) are shown in Fig. |8] where each panel cor- 
responds to different sets of values for n^ and a fixed n^. 
In each case, the individual exponents ad/on as well as 
the total exponent a are shown as a function of density. 
The exponent is extracted from a full numerical eval- 
uation of /i(n) and then using a{n) = dhi^{n)/dliin, 
where the total exponent is extracted by adding the two 
individual resistivities, i.e., /i~^ = /x^ +fj.~ ■ The rather 
amazing fact to note in Fig. [8] is that each individual ex- 
ponent ad and a^ has a maximum value ~ 1.7, albeit the 
maximum for remote (near) impurities occurring at low 
(high) carrier densities because the effective d-values are 
much higher (lower) for remote (near) impurities. (We 
note that d = for Ui impurities still has an effective d 
value of roughly a/2 sa 150 A whereas the effective d for 
the rid impurities is d -t- a/2 in each case.) Figure |8] is a 
striking direct numerical verification of our theory. 

Before discussing the same physics for graphene, which 
we do next, we mention that including screening in the 
theory is straightforward (but extremely messy). All we 
need to do is to modify Eq. (|49|) to 



/(n) 



x^e^^^^dx 



(x 4- c/y/n)'^\/l — x^ 



(55) 



with c = giving the unscreened formula of Eq. (j47l) . 
Inclusion of screening strongly affects the low-density 
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Figure 8. (color online) The calculated mobility exponent as 
a function of density for 2D n-GaAs wells of thickness a — 300 
A. The mobility is calculated by incorporating both remote 
impurity scattering with 2D impurity densities Ud and separa- 
tion d and near impurity scattering with 2D impurity density 
rii. Here we use d = 500 A and n^ = 3 x 10** cm~'^ for 
all figures, but different values of Ud, i.e., (a) Ud — 2 x 10^" 



cm"-^, (b) rid = 1 X 10^ 



(c) Ud = 0.5 X 10' 



(a) Ud = 0.22 X 10 cm~ . The calculated exponents for the 
individual scattering mechanisms a^ and a^ (for Ud and rii, 
respectively) and the total exponent a are shown as a function 
of density. 



n — >■ behavior, changing a{n — >■ 0) exponent to zero 
(for c 7^ 0) from unity (c = 0), but does not af- 
fect the intermediate or high density behavior at all (as 
is obvious from Table I where the 2k pd 3> 1 asymp- 
totic results are independent of the screening constant 
Qs ~ qTr/'^kp oc 1/y/n). Since the extrema behavior of 
interest to us is not a low-density phenomenon, our anal- 
ysis based on Eq. (j49p is appropriate (as is verified by its 
agreement with the full numerical results). 

Now, we consider the corresponding graphene case, 
which also has a maximum in a{n) at some intermedi- 
ate carrier density with Q;,nax (> 1-5) being larger than 
the corresponding infinite density unscreened exponent 
value of 3/2 (see Fig.[7|). For graphene, with the charged 
impurities located in a 2D layer a distance d from the 
graphene layer, the scattering rate t~^ is given by [see 
Eq.dall)]: 



T-i =A^/^ 



.T^yr 



dx ^\ " ^ 1 e-^"^. 



{x + qsf 



(56) 



where A, 6 = 1-^d are density-independent constants 
and qs = "ie^ /kHvo is also density independent. Direct 
expansion for Eq. ([5E)) in the low (n — >■ 0) and high {n — >■ 



-\n^0) 



Ao 



16d 

1 — VTrn 



and 



T (n — > oo) = 



Arsod 



{kpdYiqTFdY 



1- 



qrpd 



(57) 



(58) 



where Aq and A^ are constants independent of n and 
d. For graphene ^{n) ^ t / ^/n, and therefore, we get for 
n -> 



^l{n) - ( 1 + -^Vn ) , 



and for 



^JL{n) - n^/'^d^ (l + 



2rsd 



(59) 



(60) 



where rg — e'^ /(kHvo) is the so-called graphene fine struc- 
ture constant. 

Eqs. (1551) ~ (150]) imply that the mobility fj,{n) starts at 
low density (n — > 0) with a — 0, but with a leading or- 
der correction going as 0{^/n) with a positive sign. For 
large n (i.e., kpd S> 1), a{n — >■ c») becomes 3/2, but has 
a positive leading order correction of 0{l/y/n). This im- 
mediately implies that a(n) must have a local maximum 
at some intermediate density with a being only logarith- 
mically larger than the asymptotic (n — >■ cxi) value of 
3/2. This conclusion is consistent with numerical results 
presented in Fig. [71 We note that the maximum in a{n) 
for graphene is much shallower and weaker than in the 
2D parabolic system. 



VI. EXPERIMENTAL IMPLICATIONS 

We now discuss the possible experimental relevance of 
our theoretical findings in this section. The fact that 
the mobility fJ-{n) or the conductivity a{n) = nefi of 
2D carrier systems shows a density scaling behavior with 
fi ~ n°'("-) and a ~ n'^(") with P = a + 1 has been 
known for a long time in the experimental 2D transport 
literature^. Although our current work is purely theo- 
retical, focusing entirely on the fundamental questions of 
principle involving the detailed behavior of the exponent 
a{n) and /3 = a-|- 1 for various types of disorder affecting 
transport properties, we believe that it is appropriate for 
us to comment on experiments, making connection with 
the existing data in the literature as well as making some 
concrete predictions for future experiments (particularly 
in the context of our unexpected finding of a maximum 
in a{n) with an almost universal value of 1.7 for 2D semi- 
conductor quantum wells). 

We first summarize the serious difficulties in mak- 
ing direct quantitative connection or comparison be- 
tween experiment and theory in 2D semiconductor sys- 
tems with respect to low-temperature disorder-limited 
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transport properties. In fact, these caveats apply to all 
disorder-limited transport properties in all systems, not 
just to 2D transport in semiconductor quantum wells. 
The key problem is that the detailed nature of disor- 
der (either qualitative or quantitative) in a sample is 
never precisely known from independent measurements 

— in the context of our transport theory in this pa- 
per, relative amounts of 2D near and remote Coulomb 
disorder, 3D Coulomb disorder, and short-range disor- 
der are simply not known. Thus, a quantitative or even 
a qualitative theory for calculating the mobility or the 
conductivity of a given sample exists only as a matter 
of principle, but not in practice since the details of the 
underlying disorder contributing to resistive scattering 
are apriorie unknown and often are figured out indirectly 
based on quantitative comparisons between transport ex- 
periment and theory. The situation is worsened by the 
fact that the relative magnitudes of various independent 
scattering mechanisms vary strongly with carrier density 

- for example. Coulomb disorder weakens with increas- 
ing carrier density making short-range scattering rela- 
tively stronger at high carrier density. Thus, all 2D 
semiconductor transport would eventually be dominated 
by short-range scattering (e.g., interface roughness, alloy 
disorder) at very high carrier density (where a ~ in 2D 
systems) , and the only question is how high in density one 
must go to reach this asymptotic zero-range disorder lim- 
ited regime where Coulomb disorder has virtually been 
screened out. This, of course, is completely nonuniversal 
and depends entirely on the relative amount of Coulomb 
impurities and short-range disorder in particular sam- 
ples! This discussion shows that in 2D semiconductor we 
have a{n — >■ 0) = and a{n — > oo) — purely theoreti- 
cally with the zero-density limit and the infinite-density 
limit being dominated by completely screened Coulomb 
disorder {kpd ^ 1, qs » 1) and zero-range disorder, 
respectively, although these strict theoretical limits are 
unlikely to apply to real samples at any finite carrier 
density. 

The difficulty of applying the pristine theory to specific 
experimental situations is obvious from our numerical re- 
sults presented in Figs. [5] and H] (for 2D holes) and Fig.[S] 
(for 2D electrons). In each case, the mobility exponent 
a for individual near and far Coulomb impurities follows 
our theoretical prescription perfectly with the expected 
low- and high-density exponents agreeing with the results 
of Table I, but the total exponent a (which is the only one 
relevant for the experimental data) may not follow any 
well-defined pattern and could vary strongly depending 
on the relative strengths of near and far Coulomb dis- 
order, showing strong nonmonotonicity (Figs. [5] and [6]) 
or weak/no nonmonotonicity (Fig. |5|). This reinforces 
the point made earlier by us that the universal density 
scaling behavior in transport applies only to individual 
scattering mechanism with the overall transport being 
dominated by crossover behavior is generically nonuni- 
versal due to the existence of several different operational 
scattering processes. 



In spite of the above serious caveats arising from our 
ignorance about the underlying disorder contributing to 
resistive scattering mechanisms, some general statements 
can be made about the implications of our theory to ex- 
perimental data. We discuss this below. 

(i) For very dirty (and low-mobility) samples, the back- 
ground Coulomb disorder arising from the unintentional 
charged impurities in the quantum well should dominate 
transport properties, leading to a^ ^ at,. ^ 0.5 for a wide 
range of intermediate densities, (ii) When transport is 
limited by remote dopants, which would always be true 
in modulation-doped samples for kpd < 1, ae ~ 1.5 and 
ah ^ I — 1.5 depending on the hole effective mass, but 
ae/ah will decrease with increasing density as kpd > 1 
regime is reached, (iii) For modulation doped structures 
with kpd 3> 1, background disorder again dominates at 
intermediate density giving ae ^ ah ^ 0.5 — 1. 

The above situation seems to describe the existing ex- 
perimental situation for 2D quantum well transport rea- 
sonably well as discussed below. Focusing on specific 
experimental results in the literature in the context of 
our transport theory, we make the following remarks dis- 
cussing some specific experimental publications in 2D 
GaAs based electron and hole systems. 

(1) In Ref|20, the measured ae « 0.6 — 0.7 in n-GaAs 
2D system with no intentional remote dopants (the sam- 
ple is a gated undoped sample) in the density range 



10 



10 



10 cm (i.e., Qs > 1; 2kFd < 1) agrees quan- 



titatively with our theoretical results given in Figs. [31 HI 
and [8] with background 2D and 3D unintentional charged 
impurities being the main disorder mechanism as ex- 
pected for an undoped 2D system. 

(2) In a similar gated undoped 2D p-GaAs sample 
Manfra et al^ find ah ^ 0.7 for density > 10^" cm~^, 
again agreeing with our results given in Figs. [3] and U for 
background scattering. 

(3) In Harrell et al^ gated undoped n-GaAs 2D sam- 
ples, a ~ 0.6 was found for n > 10^^ cm~^ and a « 



0.33 was found for n < 5 x 10 cm 



This is both 



quantitatively and qualitatively consistent with our nu- 
merical findings in Figs. [3l 21 and [51 where scattering 
by background charged impurities in the layer leads to 
a w 0.3 — 0.7 in the n — 10^" — 10^^ cm~^ density range 
with a{n) decreasing with decreasing carrier density. 

(4) Melloch et alM. found a « 0.6 - 0.7 for n > 10" 
cm~^ which is consistent with our background impurity 
scattering results. 

(5) Pfeiffer et al. studied^! modulation-doped high- 
mobility 2D GaAs electron systems^'* obtaining a ^ 0.7 
around n ~ 3 x 10^^ cm~^ for modulation-doped struc- 
tures {d = 1000 - 2000A) with /i > 10^ cmVVs. Again, 
remote impurity scattering is completely ineffective here 
because kpd ^ 1 rendering mobility limited only by re- 
mote impurity scattering to be around 10^ cm^/Vs ac- 
cording to our numerical calculations. The dominant 
scattering mechanism in this sample is by background un- 
intentional charged impurities, leading to a w 0.7 around 
n '-^ 3 X 10^^ cm~^ according to oiu Fig. [H which is in 
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precise agreement with the date of Pfeiffer et al,^'^. 

(6) In a similar high-mobility modulation-doped 2D n- 
GaAs sample, Shayegan et al^ found a w 0.6 in samples 
with n « 10^ cm^/Vs for n < 10^^ cm~^ with the spacer 
thickness d = 1000 — 2000 A. Again, remote scattering by 
the intentional dopants is ineffective as a resistive scat- 
tering mechanism here with the dominant scattering be- 
ing by unintentional background impurities in the GaAs 
quantum well. From Fig. |4] of our presented results, we 
find a « 0.6 for n < 10^^ cm~^ in agreement with the 
experimental finding of Shayegan et al^. 

(7) Most of the high-mobility experimental samples 
discussed above are dominated by the background unin- 
tentional charged impurities in the 2D layer itself leading 
to a < 1 by virtue of the fact that the intentional dopants 
introduced for modulation doping are rather far away in 
these high quality samples (this is a generic feature of all 
high mobility 2D samples with ji > 10^ cm^/Vs where 
a < 1 prevails by virtue of the background disorder be- 
ing dominant). By contrast, early work on modulation- 
doped 2D samples invariably had lower values of d and 
achieved much lower mobility /u < 10® cm^/Vs. Such 
samples are almost always dominated by remote scatter- 
ing by the intentionally introduced dopants, leading to a 
values typically exceeding unity as our theory predicts. 
As a typical example, we consider the work of Hirakawa 
and Sakaki2£ on modulation-doped 2D n-GaAs samples 
with /i - 10" - 5 X 10^ cmVVs for d « - 180 A in 



the 



10^ 



5 X 10 cm density range. Assuming 



transport to be limited entirely by the intentional ionized 
dopants in the modulation layer (i.e., no background dis- 
order scattering), our results of Fig. [T] predict awl for 
d = and a « 1.1 — 1.3 for d = 100 A. Hirakawa and 
Sakaki reported^^ a w 1.1 - 1.3 for d w - 100 A in 
essential agreement with our theory. In addition, these 
authors reported a « 1.7 for d « 200 A, an anomalous 
mobility exponent (i.e. d > 3/2) which has remained un- 
expected in the literature for more than 25 years. Our 
current work provide a definitive explanation for a w 1.7 
as arising from the maximum in a{n) for remote scatter- 
ing as is apparent in Fig. [1] We note that kpd '^ 1 — 2 for 
n ~ lO'^^ cm~^ and d ^ 100 — 200 A, and thus our theory 
predicts a ~ 1.7 for d sa 200 A in the Hirakawa- Sakaki 
experiment^®. We believe that the experimental finding 
of a ^ 1-7 by Hirakawa and Sakaki is a direct verification 
of our intriguing prediction of a > 3/2 for kpd > 1 in 
transport dominated by remote scattering. 

(8) Finally, we discuss some recent unpublished exper- 
imental work by Pfeiffer and West^^ who, motivated by 
our theoretical work, carried out low-temperature trans- 
port measurements in a series of high-quality (i.e. ultra- 
pure GaAs with very little background disorder due to 
unintentional impurities) MBE-grown 2D n-GaAs sam- 
ples with variable values of d. Since these experiments 
were performed with the specific goal of checking our 
low-temperature 2D transport theory predictions, Pfeif- 
fer and West made undoped gated samples of highest 
quality with little background disorder and a nominal 



low-temperature mobility oi fi > lO"^ cm^/Vs. Then, 
they systematically introduced charged impurities at spe- 
cific separation (d) from the 2D layer by inserting carbon 
atoms in the GaAs layer (d = 0) or in AlGaAs barrier 
layer (d ^ 0). First they explicitly verified that the in- 
troduction of different amounts of impurity centers with- 
out changing d only affects the 2D mobility through the 
expected rii^ scaling behavior (i.e. fi ex n~ ) without 
changing a{n) by changing the inserted carbon atoms 
by a factor 5 keeping d fixed (which changed the 2D 
mobility by a factor of 5 without changing the expo- 
nent a{n) in the same carrier density range). Thus, 
they measured a{n) for d — and d = 150 A finding 
a{n) « 0.8 and 1.8 respectively for n ~ 10^^ cm^^. Our 
calculated a{n « 10" cm-^)^ 0.8 for d = in Fig. [2] 
in perfect agreement with the experimental data. For 
d — 150 A, kpd sa 1 for n w 10^"'^ cm~^, and we predict 
a{n) sa a^a.^ = 1.7 in this situation which compares well 
with the experimental exponent of 1.8. Thus, this exper- 
imental investigation of our theoretical predictions ap- 
pears to have strikingly verified our finding that a > 3/2 
in the kpd ^ 1 intermediate density regime where the 
shallow maximum occurs in a{n) in our theory. 

(9) Before concluding our discussion of experimental 
implications of our theory we describe the very recent 2D 
hole transport data in high-mobility p-GaAs systems by 
the Manfra group ^^'^^ . These 2D p-GaAs samples with 
hole mobility > 2 x 10® cm^/Vs are the world's high- 
est mobility hole samples ever, and taking into account 
the effective mass difference (~ a factor of 5) between 
GaAs electrons and holes compare favorably with the 
best (^ 15x 10® cm^/Vs) available 2D electron mobilities. 
The main finding of the work^ - is that the mobility expo- 
nent a{n) increases from 0.7 at high hole density (^ 10^^ 
cm^2) to 1.7 at low density (> 10^" cm~^) for modula- 
tion doped sample with d — 800 A. Since the mobility 
remains high throughout (^ 10^ cm^/Vs), localization 
effects should not be playing a role. We therefore believe 
that the experimental finding of Watson et ali^ is a direct 
confirmation of the a{n) behavior for 2D holes presented 
in our Figs. [5] and El where the total calculated a{n) in- 
creases monotonically as the 2D hole density decreases 
from n ^ 10^^ cm~^ to n ^ 10^^° cm~^. In fact, even the 
explicit a-values measured by Watson et al. agree well 
with our 2D hole theoretical results presented in Figs. [5] 
and [6] with a (n ~ 10^*^ cm~^) increasing above the un- 
screened a = 3/2 value reaching essentially the measured 
value of a ^ 1.7 around n > 10^" cm^^. Thus, our the- 
ory provides a qualitative explanation of the Watson et 
al. experimental results including the surprising finding 
of the intermediate-density a being around 1.7 (> 1.5). 

(10) We conclude this section on the experimental 
relevance of our theory by discussing graphene briefiy. 
There has been substantial research activity on studying 
the density dependent graphene conductivity^.— which 
is well beyond the scope of our current work and has al- 
ready been covered elsewhere. It is known that scattering 
by near random charged impurities located on the surface 
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of the graphene layer or at the graphene-substrate inter- 
face leads to a (7{n) ex n (i.e., fi{n) ^ constant) in the 
intermediate density {kpd < 1) region, and in the high- 
density regime a{n) becomes sublinear most likely be- 
cause of short-range defect scattering which gives (t{ti) ^ 
constant (i.e., ^ ^ ^/n) - see Table I for details. Theory 
predicts (Table I) that for kpd ^ 1, i.e., remote Coulomb 
disorder, /i(n) and a{n) should cross-over to /i(n) ~ ^/n 



and cr{n) 



7.3/2 



in graphene. This clear prediction could 



be verified by putting an impurity layer (e.g. a Si02 
film) at various values of d from the graphene layer and 
measuring a{n) to check if the low-density (2fci?d ^ 1) 
a{n) ^ n behavior indeed crosses over to the high den- 
sity (j{n) ^ n^" behavior as we predict theoretically. In 
graphene, with a valley degeneracy of 2, 2fci?(i ss 2.5d\/n, 
where d — d/lOOOA, and h — n/10^°cm^^. Hence, for 
n = 10^2 cm-2 and d = 100 A 2kFd « 2.5. Thus, 
the condition 2kpd 3> 1 would require an impurity layer 
at d sa 1000 A (with consequent very weak Coulomb 
disorder scattering) which may lead to the complication 
that the subsequent graphene resistivity will be entirely 
dominated by any underlying short-range disorder (with 
a ^ n^), masking any a ^ rc'l'^ behavior arising from 
charged impurity scattering. One possibility would be to 
put suspended graphene near a thick layer (of thickness 
L) of disordered substrate with a 3D charged impurity 
distribution, which would lead to [see Eqs. (|38p ~ pn)) : 



A^, 



nh^vok^ 



2TTe' 



dx- 



;%/! 



[x + qsY 



[1 



-2LoxT, 



(61) 

with Ni being the 3D impurities density in the impurity 
layer and Lq = '^kpL, giving cr(ri) ~ n^'^ for kpL ^ 1 
and cr{n) ^ n for kpL <C 1. Such a 3D impurity 
layer underneath graphene may manifest the superlinear 
(j{n) ~ n^'"^ conductivity behavior predicted by the the- 
ory, but observing the shallow maximum in the exponent 
PI a for kpd ^ 1 may still be difficult. 



VII. OTHER EFFECTS 

In this section, just before our conclusion in the next 
section, we discuss "other effects" completely left out 
of our theoretical considerations which may compromise 
and complicate direct quantitative comparisons between 
our theory and experiment although we believe that our 
theoretical conclusions should apply generically to trans- 
port in high-mobility 2D semiconductor systems at low 
enough temperatures. 

First, phonon effects are neglected in the theory since 
we explicitly consider the T = situation (in practice, 
T = 50 — 300 mK is the typical low-temperature exper- 
imental situation mimicking the T = theoretical situ- 
ation). For consistency between theory and experiment, 
the transport data must therefore be taken at a fixed 
temperature T < Tbg, where Tbg is the so-called Bloch- 
Griineisen temperature, so that acoustic phonon scat- 



tering contribution to the resistivity is negligible com- 
pared with disorder scattering even in the high-mobility 
2D semiconductor structures under consideration in this 
work. (Optical phonon transport is of no relevance for 
low-temperature transport since ksT <^ fuj->LO is explic- 
itly satisfied as huj^o > 100 K typically.) Tbg is given 
either by the Debye temperature (for 3D metals) or by 
the energy of the acoustic phonons with 2fci?-wave vector 
(for 2D semiconductors), whichever is lower. We there- 
fore have kBTBG — '2'fikpVph where Vph is the relevant 
acoustic phonon velocity in the material. Putting in the 
appropriate sound velocity {Vph), we get Tbg ~ 2vn K; 
10\/n K in 2D GaAs and graphene, respectively, where 
fi is measured in units of 10^° cm~^. Thus, down to 
carrier density n ^ 10^ cm~^, it is reasonable to ignore 
phonon effects in transport at T w 100 mK. Acoustic 
phonon scattering has been considered elsewhere in the 
literature^. 

Second, we have used the Born approximation in cal- 
culating the scattering time and the RPA in calculat- 
ing the screened Coulomb disorder throughout. Both of 
these approximations surely become increasingly quanti- 
tatively inaccurate at lower carrier density although they 
should remain qualitatively valid unless there is a metal- 
insulator transition (obviously, our Drude-Boltzmann 
transport theory would not apply at or below any metal- 
insulator transition density). RPA screening theory be- 
comes increasingly quantitatively inaccurate as carrier 
density (or the corresponding r^ ~ n~^'^ parameter) de- 
creases (increases), but there is no well-accepted system- 
atic method for incorporating low-density electronic cor- 
relation effects going beyond RPA screening. Including 
low-density correlation effects in Hubbard-type local field 
theories do not change any of our qualitative conclusion. 
As for multiple scattering effects™ beyond Born approx- 
imation, they typically lead to higher-order corrections 
to the resistivity so that p on nt, where nt is the impu- 
rity density, is no longer valid and one must incorporate 
high-order nonlinear corrections to the resistivity going 
as 0{nf) and higher. These nonlinear multiscattering 
corrections become quantitatively important for n "^rii, 
and can be neglected for n > Ui regime, which is of our 
main interest in this paper. Our neglect of multiscat- 
tering corrections beyond Born approximation is consis- 
tent with our neglect of strong localization effect, both of 
which will become important in the very low carrier den- 
sity regime {n < rii) where the Drude-Boltzmann theory 
becomes manifestly inapplicable. 

Third, we ignore all nonlinear screening effects, which 
have been much discussed in the recent graphene 
literature^ where charged impurity induced inhomoge- 
neous electron-hole puddle formation plays an impor- 
tant role at low carrier density, since they are impor- 
tant only at very low carrier density (n < rii) where our 
whole Boltzmann theoretical approach becomes suspect 
any way. For the same reason we also do not take into 
account any scattering-screening self-consistency effeclj^ 
which may also become important at very low carrier 
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density (again, for n < rii). 

Fourth, we ignore any possible spatial correlation ef- 
fects among impurity locations, assuming the disorder to 
arise from completely uncorrelated random impurity con- 
figurations. If the impurities are spatially correlated, it 
is straightforward to include the correlation effect in the 
Boltzmann transport calculation by simply multiplying 
the disorder potential term in Eq. (|5]) by the correspond- 
ing structure factor s{q) for the impurity distribution, 
i.e., by writing \Vq\'^ — ^ |T/jps(g) in Eq. ^ where 



n,; 



E' 



n,5, 



iOqO, 



(62) 



where ri denotes the position of each impurity and the 
sum going over all the impurities. This, if experimen- 
tal information about impurity correlations exists, it is 
then straightforward to include such spatial correlation 
effects in the Boltzmann transport calculations, as has 
indeed been done in both 2D GaAs"^"* and graphene'^-. 
Since intrinsic impurity correlation information is typi- 
cally unavailable, our model of uncorrelated random dis- 
order seems to be the obvious choice from a theoretical 
perspective since it involves only one (the impurity den- 
sity rii) or two {rii and the impurity location d) unknown 
parameters (and most often d is known from the mod- 
ulation doping setback distance and is thus not a free 
parameter) whereas including impurity spatial correla- 
tions would invariably involve the introduction of more 
unknown free parameters making the theory of dubious 
theoretical relevance. Since our interest in this paper is 
not an absolute calculation of the conductivity or mo- 
bility (and hence ni typically drops out of our theory if 
only one dominant scattering mechanism is operational) , 
but obtaining the universal density dependence of con- 
ductivity, it is important to mention that the main effect 
of impurity correlation is to suppress the effect of n,; on 
the resistivity without much affecting the carrier density 
dependence particularly for n > Ui regime of our interest. 
We mention that any intrinsic impurity correlations actu- 
ally increase the value of mobility from our uncorrelated 
random disorder theory and thus compensate to some 
extent the suppression of mobility arising from some of 
the effects discussed above. 

Fifth, for our T — theory to be strictly applicable to 
the experimental data, the experimental temperature T 
should satisfy the condition T <C Tp, where Tp = Ep/ks 
is the corresponding Fermi temperature of the system. 
Using the known dependence of Ep on the 2D carrier 
density n we see that this implies T{K) <^ A.2n for 
2D n-GaAs, T{K) < In for 2D p-GaAs, and T{K) < 
150\/n for graphene where fi = n/10^'^cm~^. Thus, low- 
temperature experiments carried out at T = 100 mK 
satisfy our T — theoretical constraint very well down 
to 10^ cm~^ carrier density (except for 2D p-GaAs hole 
system where the density cut off is perhaps 5x 10^ cm~^). 
Our direct numerical calculations at finite temperatures 
(not shown in this paper) show that our calculated ex- 



ponents a(n) and /?(n) =a + latT=:0 continue to 
be quantitatively accurate upto T > Tp as long as the 
experimental data for density dependence are taken at 
a fixed temperature. Thus our theory and numerics for 
a{n) and /3{n) are quite robust against thermal effects 
as long as phonons are unimportant (i.e., T < Tbg is 
satisfied). 

We conclude this section of "other effects" by noting 
that the most important drawback of our RPA-Drude- 
Boltzmann theory is that it may fail systematically at low 
carrier density (n < rii) where important physical effects 
(which are difficult to treat theoretically) such as strong 
localization, metal-insulator transition, nonlinear screen- 
ing, inhomogeneous puddle formation, multiscattering 
corrections, screening-scattering self-consistent, etc. may 
all come into play making both our theory inapplica- 
ble to the experimental situation. We do, however, an- 
ticipate our theory to be applicable to very low car- 
rier densities (n > 10^ cm~^) in ultra-high mobility 2D 
GaAs and graphene systems where rii ^ 10^ cm^^ typ- 
ically. A convenient experimental measure of the ap- 
plicability of our theory is looking at the dimension- 
less quantity ^kpV where I is the elastic mean free path 
defined by I = vpT where r and vp are respectively 
the transport relaxation time and Fermi velocity. We 
find that kpl = 4.14n/i for 2D GaAs systems, where 
fi = n/10^°cm^^ and jl = /j,/(10® cm^/Vs) whereas 
kpl ~ 0.2n/i for graphene. As long as kpl > 1, our 
Drude-Boltzmann theory should be valid qualitatively 
and we therefore conclude that the theory remains quan- 
titatively accurate for n > 10^" cm~^ (10^^ cm~^) for 2D 
GaAs (graphene) systems in high-quality/high-mobility 
samples. Thus, there is a large range of carrier density 
(10^° - 10^2 for 2D GaAs and 10" - 10" for graphene) 
where our predictions for universal density dependence 
can be experimentally tested through low-temperature 
transport measurements. 

The fact that all the "other effects" left out of our 
theory affects transport at low carrier densities indicates 
that our predicted density scaling of conductivity will 
systematically disagree with the experimental data at 
lower carrier densities. Of course, "lower" density is a 
relative term, and the dimensionless quantities such as 
n/rii and kpl are the appropriate quantities to define the 
regime of validity of our theory. As n/rii and/or kpl 
(or even Tp/T or Tbg/T) become smaller, the Drude- 
Boltzmann theoretic predictions become increasingly un- 
reliable. Nevertheless, the theory remains predictive 
down to 10^*^ cm~^ carrier density (or lower) in high 
quality (i.e. low values of rii) GaAs samples at low tem- 
peratures (Tsa 100 mK). 



VIII. DISCUSSION AND CONCLUSION 

We have developed a detailed quantitative theory for 
the density-dependence of the zero-temperature conduc- 
tivity (or equivalently mobility) of (mainly) 2D and 3D 
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electron and hole metallic systems assuming transport 
to be limited by (mainly Coulomb) disorder scattering 
within the semiclassical Drude-Boltzmann transport the- 
ory. We neglect all quantum interference (hence localiza- 
tion) effects as well as interaction effects (except for the 
carrier screening of the bare impurity Coulomb disorder, 
which is an essential qualitative and quantitative ingre- 
dient of our theory) assuming them to be small since 
our interest is the density dependence (rather than the 
temperature dependence) of transport properties at low 
fixed temperatures in high-mobility {kpl 3> 1 where I is 
the elastic mean free path) samples. 

We have systematically considered 3D and 2D doped 
(n- and p-) semiconductor systems as well as 2D graphene 
but the primary focus has been on n-GaAs and p-GaAs 
quantum well based 2D electron or hole systems, mainly 
because these systems continue to be of great interest 
in physics and because the carrier density can easily 
be tuned in such high-mobility 2D semiconductor sys- 
tems, and the mobility is dominated by Coulomb disor- 
der at low temperatures. We have taken into account 
both long-range Coulomb disorder from charged impu- 
rities and zero-range disorder arising from possible non- 
Coulombic short range scatterers. The primary focus has 
been the Coulombic disorder since this is the main low- 
temperature resistive scattering mechanism in semicon- 
ductors. Instead of discussing the nonuniversal values of 
/U (and a), which depend^^ on the actual impurity con- 
tent, we focus on the universal power-law density scaling 
of transport properties: cr(n) ^ n^M and fJ.{n) ^ n"^"' 
with /3 = a + 1. These exponents a and /? are sample- 
independent and depend only on the nature of the domi- 
nant disorder. We provide asymptotic theoretical analy- 
sis of a and /? (for various types of underlying disorder) in 
the high and the low density regime and for near and far 
impurities. We have then verified our analytical results 
with direct numerical calculations based on the full solu- 
tion of the Boltzmann transport theory in the presence 
of disorder scattering. 

Although our work is primarily theoretical, we pro- 
vide a critical comparison with various experimental re- 
sults in the literature (in 2D n-GaAs electrons and p- 
GaAs holes), finding generally good agreement between 
our theoretically predicted exponents a and /3 and low- 
temperature experimental findings in high-mobility 2D 
electron and hole systems. In particular, a truly excit- 
ing prediction, that a{n) has a maximum universal value 
cimax ^ 1-7 for all 2D systems at an intermediate car- 
rier density value (approximately defined by kpd ^ 1), 
seems to be consistent with recent (and old) experimen- 
tal results from several different groups (and for both 2D 
electrons and holes), as discussed in details in Sec. VI. 

Our theory predicts a{n) to vary from a = in the 
low-density strong-screening regime (qte ^ 2fci?, n — > 0) 
to a = 3/2 in the high density weak screening regime 
ilTF <^ 2fcF, n — > oo) with a shallow maximum of 
a ~ 1.7 at intermediate carrier density for kpd ^ 1 
where d is the impurity location. Although our work pre- 



sented in this article is purely theoretical describing the 
density dependent and disorder-limited T — Q conductiv- 
ity of 2D/3D carriers using the semiclassical Boltzmann 
theory approach, it is worthwhile to speculate about the 
prospects for the experimental observation of our asymp- 
totic low (a ->; ^ -> 1) and high [a -^ 3/2; /3 -^ 5/2) 
density behavior of Coulomb disorder limited 2D semi- 
conductor transport. These limiting exponents are theo- 
retically universal. 

We first summarize the current experimental status in 
the context of our theory. For scattering by background 
Coulomb disorder (i.e., near impurities with 2kpd < 1), 
a{n) S3 0.5 — 0.8 typically, and for scattering by remote 
impurities (2kpd > 1) in the modulation doping layer 
a{n) > 1 — 1.3 typically with a few atypical cases show- 
ing a{n) ~ 1.7 (> 3/2) around kpd ~ 1. All of these 
are intermediate density results in our theory where qs 
{= qTp/2kp) and 2k pd are neither extremely large nor 
extremely small. Thus the basic experimental situation 
is in excellent agreement with our theory as it should 
be because the 2D doped semiconductor transport (as 
well as graphene) is known to be dominated by screened 
Coulomb disorder with near or far charged impurities be- 
ing the dominant scattering mechanism depending on the 
sample and carrier density. 

First, we discuss the high-density situation which is 
theoretically more straightforward and where the Boltz- 
mann theory is almost exact. As carrier density increases, 
the semiclassical Boltzmann theory becomes increasingly 
more valid for Coulomb disorder limited transport prop- 
erties since the conductivity itself and consequently kpl 
increases, thus making the system progressively more 
metallic. In addition, increasing density decreases the 
metallic r^ -parameters, the dimcnsionless Wigner-Seitz 
radius, given by r^ = rae^ j {kTi} \pKn) for 2D semicon- 
ductor systems. Since r^ ^ n~^''^ , at high carrier density 
(e.g., rg = 0.5, 2.5 for 2D n-GaAs, p-GaAs respectively at 
n = 10^^ cm~^) Tg is small, making our theory using RPA 
screened effective Coulomb disorder to be systematically 
more valid at higher carrier density as RPA becomes ex- 
act at low Ts- Thus, it appears that the ideal applicability 
of our theory is in obtaining the high-density 2D system. 

This is indeed true except that new physical (rather 
than theoretical) complications arise making it problem- 
atic for a direct comparison between our theory and ex- 
periment on 2D systems at high carrier density. Two new 
elements of physics come into play at high carrier den- 
sity, both contributing to the suppression (enhancement) 
of mobility (scattering rate): Intersubband scattering be- 
comes important as the Fermi level moves up and comes 
close to (or crosses over into) the higher confined sub- 
bands of the quasi-2D quantum well, thus opening up 
a new scattering channel; and, short-range scattering at 
the interface and by alloy disorder (in AlGaAs) becomes 
important as the self-consistent electric field created by 
the electrons themselves pushes the carriers close to the 
interface at high carrier density. Both of these physi- 
cal effects eventually suppress the monotonic growth of 
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li{n) and a{n) with increasing density, and eventually 
fi{n) starts decreasing with increasing carrier density at 
high enough density (for n > 3 x 10^^ cni^'^ in GaAs- 
AlGaAs systems) instead of continuing as ii{n) ^ r?^"^ 
as it would in the high-density regime if Coulomb disor- 
der is the only dominant scattering mechanism with no 
other complications. Since at high density r^^(n) ^ n" 
with a > 1 for Coulomb disorder, eventually at some 
(nonuniversal sample dependent) high density. Coulomb 
disorder becomes insignificant compared with the short- 
range scattering effects. Obviously, the density at which 
this happens is non-universal and depends on all the de- 
tails of each sample. But in all 2D systems and sam- 
ples, eventually, when Coulomb disorder induced scat- 
tering rate is sufhciently small, a high density regime is 
reached where the mobility stops increasing (and even 
starts decreasing). In Si MOSFETsi^ this effect is very 
strong already around ^ 10^^ cm~^ because of consid- 
erable surface roughness scattering at the Si-Si02 inter- 
face, and [i{n) decreases with increasing n at higher den- 
sity. Even in high-mobility GaAs systems, /i(n.) saturates 
(i.e., a = 0) and eventually starts decreasing at a non- 
universal density around J^ ^ 3 x 10^^ cm~^. Boltzmann 
transport theory can be easily generalized to incorporate 
inter-subband scattering and surface scattering, but the 
physics is nonuniversal and beyond the scope of our cur- 
rent work. 

The low-density (n — >■ 0) situation is fundamentally 
inaccessible to our Boltzmann transport theory since all 
doped semiconductor systems (3D or 2D) eventually un- 
dergo an effective metal-insulator transition at (a nonuni- 
versal) low (critical) density with the semiclassical Boltz- 
mann theory eventually becoming invalid as a matter 
of principle at a sufficient low sample- dependent carrier 
density. In 3D, this transition may be a true Anderson lo- 
calization transition as kpl ^ ^ with decreasing density 
making the Boltzmann theory inapplicable. In 2D semi- 
conductors, which are the systems of our main interest, 
the observed metal-insulator transition at a nonuniver- 
sal sample dependent critical density ric is likely to be a 
crossover phenomenon^'^ since the scaling theory of An- 
derson localization predicts 2D to be the critical dimen- 
sion with no localization transition. There is considerable 
experimental support for the observed low-density 2D 
metal-insulator transition to be a density inhomogeneity- 
driven percolation transition at a nonuniversal critical 
density Uc where charged impurity induced Coulomb dis- 
order drives the system into an inhomogeneous collection 
of "puddles" with a mountain and lake potential land- 
scape where semiclassical metallic transport becomes im- 
possible for n < Tic with Uc being the disorder-dependent 
percolation transition densityS2ii^i^"— . The critical den- 
sity Tic depends crucially on the sample quality and is typ- 
ically below 10^ cm^^ in high-mobility GaAs systems — 
in Si MOSFETs, where disorder is very strong, ric « 10^^ 
cm~^, and this is why we have left out 2D Si systems from 
our consideration in this work although our basic theory 
applies well to 2D Si-based systems for n > 10^^ cm~^. 



Our semiclassical Boltzmann theory works for n ^ Uc, 
but for n ^ Uc, one must include the inhomogeneous 
puddle formation and the associated percolation transi- 
tion even in the semiclassical theory"^^. This is the reason 
our theory fails for real 2D systems in the n — >■ limit 
unless the level of disorder is extremely low (even then 
our Boltzmann theory is valid only for n ^ ric and the 
n — )• limit is fundamentally inaccessible). 

In spite of this fundamental difficulty in accessing the 
low-density (i.e. n — >■ 0) limit using our theory directly, it 
turns out that we can approximately include the effect of 
a semiclassical percolation transition in our theory by an 
indirect technique in the density regime n > ric above the 
effective metal-insulator transition (but still in a reason- 
ably low density regime for high-quality samples where Uc 
can be very low) . Since the percolation picture essentially 
eliminates a certain fraction of the carriers from being 
metallic, we can assume that the effective conductivity 
(mobility) for n 3> nc is given by the same exponent /3 
(a) calculated in our theory with the only caveat that /3 
{a) is now the exponent only for the actual "metallic" 
free mobile carrier fraction of the whole system. We can 
then write: 



/i = A{n - Uc)" = Bn" , 



(63) 



where a(n) is the real exponent we calculate theoretically 
from the Boltzmann theory and a'{n) is the effective (ap- 
parent) exponent obtained experimentally from 



a In) 



dlnfi{n) 
dlnn 



(64) 



by using the actual data for fi{n) without taking into 
account any complications arising from the existence of 
Tie- Using Eqs. (|551) and (IMl) . we immediately get the 
following relationship connecting the effective exponent 
a' with the real exponent a for n'^ ric 



(65) 



a(I - Uc/nY 



for n^ ric. We note that a' « 2a for n = 2?ic- Eq. (|65|) 
is valid within logarithmic accuracy, and connects the 
measured low-density {n small, but n ^ iic being still 
valid) exponent a' [n) with the real Boltzmann theory 
exponent (as in Table I) a{n). Since ric can be mea- 
sured experimentally by checking where (j{n) vanishes, 
and defining ric to be a(ric) = at T = 0, we then im- 
mediately see that the measured low-density effective ex- 
ponent a'(n) > a{n) always, and in fact, a'{n) -^ oo 
as n — >■ ri^ . Thus, we conclude that the experimen- 
tally measured mobility exponent will eventually start 
increasing as n decreases approaching ric- For n '^ ric, 
we get a' M a, but the leading correction to a goes as 
a' = a(I + nc/n) for n ^ ric- We note that the existence 
of ric enhances a over its nominal value of Table I even 
for n 3> ricl 

We have checked existing experimental results in the 
literature (to the extent ric, ct' etc. are known experi- 
mentally), finding that our prediction of Eq. (1651) seems 
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to apply quite well. We note that since the pristine cal- 
culated a{n) decreases with decreasing n (see Figs. [T]-[5]) 
for all models of Coulomb disorder at low carrier density, 
Eq. (1651) predicts that the apparent exponent a' would 
be close to a for n ^ ric, but would then manifest a min- 
imum at a density no (> ric) defined by the equation: 



da a{nc/n) 
dn {n — fic) 



(66) 



at n = no > n^,. For n < ng, a' will increase as n 
decreases (with a' being a minimum at n = no), even- 
tually diverging as (1 — nc/n)~^ as n -^ n+. We find 
this behavior to be qualitatively consistent with all the 
existing data in the literature although a precise quan- 
titative comparison necessitates more low temperature 
data showing /^(n) all the way down to n — ric where 
fi{n) = 0. Careful measurements of cr(n) close to Uc are 
lacking in the literature for us to form a definitive con- 
clusion on this matter at this stage. 

We conclude by pointing out that very close to the per- 
colation transition, where n — nc <C n^, i.e. n <C 2nc, we 
expect the critical behavior of 2D percolation transition 
to possibly come into play, where the conductivity may 
have a completely different universal 2D percolation ex- 
ponent S (totally distinct from a or a') which has nothing 
to do with our Boltzmann theory: 



(n - n,f 



(67) 



for (n — Uc) <^ ric- This percolation critical exponent 
5 for cr(n) can only manifest itself very close to n^ (i.e. 
n ^ 2nc), and for n > 2nc we believe that our effective 
theory predicts an effective conductivity exponent /3' ~ 
a' -\-l given by 



P' 



1 — Tic/n 



1 



P - nc/n 
1 — ric/n ' 



(68) 



for n > 2nc. Again, as for a', /?' w /3 = a + 1 for n ^ Uc- 
The conductivity exponents S and /3' arise from com- 
pletely different physics (from percolation critical theory 
near ric and Boltzmann theory far above Uc^ respectively) 
and have nothing to do with each other. How the tran- 
sition or crossover occurs even within the semiclassical 
theory from (3' (for n > 2nc) to 6 for (n ^ 2nc) is a very 
interesting question which is beyond the scope of the cur- 
rent work. In our current work, we have, however, solved 
[Eqs. (|63|) - (|68)) above] the problem of the crossover from 
the effective exponent a' (or (3') for n > 2nc to the true 
Boltzmann exponent a (or /3) for n 3> 2nc. 

We give one possible experimental example for the ob- 
served crossover from the Boltzmann exponent a (for 
n 3> ric) to a' (for n > ric) arising in Jiang et al^. In this 
work^, the conductivity was measured in a high-mobility 
modulation doped 2D GaAs electron system, finding the 
effective mobility exponent a « 0.9— 1.1 in the high den- 
sity range (n - 5x lO^^-Sx 10" cm'^) for d = 350-750 
A. This range of n and d converts to 2kFd =1 — 5 and 



(7s ~ 1, which is the intermediate density range for all 
our Coulomb disorder mechanisms in Table I. Given that 
the measured mobility in Ref. [40] was relatively mod- 
est, /i ^ 10^ — 10^ cm^/Vs, it is reasonable to expect, 
based on our numerical mobility calculations, for both 
remote and background Coulomb scattering to be equiv- 
alently effective, leading to a ~ 0.9—1.2 according to our 
numerical results of Sec. IV. Thus, the "high-density" ex- 
ponent (a ~ 1) in the Jiang sample agrees well with our 
theory. The interesting point to note in the current con- 
text is that Jiang et al^ found a large increase of the 
measured exponent a as n decreases, which is in appar- 
ent disagreement with our analytical theory which always 
gives a{n) decreasing with decreasing n. Although other 
possibilities (e.g., multiscattering) cannot be ruled out^i, 
we believe that the Jiang et al. data indicate a crossover 
from a to a' as n decreases. In particular, the mea- 
sured mobility exponent increased from a ^ 1 to a ^ 4 



as n decreased from 



to 



2 X 10^^ cm' 



(with also a factor of 50 decrease in mobility). Assum- 
ing ric ~ 1-5 X 10^^ cm~^ (which is consistent with the 
data), we get [from Eq. ([53])], a' w 4a Ri 4 at n « 2 x 10^° 
cm~^, which is in excellent agreement with the data^. 
We therefore believe that the Jiang et al. experiment 
manifests our predicted crossover behavior from a{n) to 
a'{n) as n approaches ric from above. 

We note that kpl ^ 1 in the Jiang et al. sample for 
n w 3 X 10^° cm~^ (where fi sa 10^ cm^/Vs seems to 
have been reached starting from /i ~ 10^ cm^/Vs for 
n ~ 3 X 10^^ cm~^), and thus the identification of ric ~ 
1.5 X 10^" cm^^ is meaningful. Clearly, our simple Drude- 
Boltzmann theory applies for n ^ 3 x lO^*' cm~^, but 
not for n < 3 X 10^" cm~^ where kpl ^ 1. What is 
encouraging is that the simple modification of the theory 
introducing the crossover exponent a' ~ a{l — n/nc)~^ 
seems to describe the experimentally observed density 
scaling exponent of the observed experimental mobility 
at low carrier densities. 

We conclude by mentioning that direct numerical per- 
colation calculations indicate 6 ^ 1.32 in 2D systems, 
which is unfortunately too close to our calculated value of 
/? in the low-density Coulomb disorder-dominated Boltz- 
mann theory, where (see Figs.[l]-|6]and[8l) the low-density 
(n < 10^*^ cm~^) a-value is a « 0.3 — 0.5 implying low- 
density /3 « 1.3 — 1.5. Since high-mobility 2D GaAs sam- 
ples typically have nc « 10^ cm~^, it is unclear whether 
the existing conductivity exponent measurements for the 
putative metal-insulator transition2£ii^i^i22 for the 2D 
density-driven metal-insulator transition really is obtain- 
ing S or is just a measurement of our calculated /3 = a-\-l 
which at low values of n would be rather close to the ex- 
perimentally measured percolation exponent of 1.3— 1.5. 
The current experiments do not really quantitatively 
measure (T(n) in the n < 2nc regime necessary for obtain- 
ing S, and we feel that much more work will be needed to 
establish the nature of the a{n — >■ ric) ^ transition ob- 
served in the laboratory. It is possible, even likely, that 
the existing measurements have only measured the low 
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density (but still n ^ ric) value of our calculated (non- 
critical) Boltzmann exponent /? « 1.3 — 1.5 in n- and p- 
2D GaAs systems. 



theoretically and numerically, the power law exponents 
for the density-dependent mobility and conductivity 
within the Boltzmann transport theory. 



In conclusion, we have developed a comprehensive 
theory for the universal density scaling of the low- 
temperature transport properties of 2D and 3D doped 
semiconductors and graphene, concentrating on the role 
of background Coulomb disorder and obtaining, both 
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